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Abstract. In order to determine the range of validity of certain elementary concepts 
in waveguide theory, the propagation of electromagnetic waves along and through a 
coaxial line with imperfectly conducting walls is studied in some detail for a particular 
method of driving. In particular, it is found that the usual concept of attenuation is 
meaningful only for a certain range of distances from the driving point. Beyond this 
distance, the electromagnetic field in the coaxial line behaves more like a radiation 
field. The explanation is supported by the behavior of the electromagnetic field in the 
imperfect outer conductor of the coaxial line. It is also found that the solution in terms 
of the ‘“‘mode” concept has a surprisingly limited region of validity. The reflection 
coefficient, the radiation pattern and the transverse distribution are also determined. 

1. Introduction. The propagation of a wave through an imperfectly conducting 
coaxial line is of interest from the point of view of waveguide theory. In the elementary 
theory of waveguides, the wall of the waveguide is assumed to be perfectly conducting. 
Under this assumption, there arises the extremely important, useful, and convenient 
concept of “mode.’”’ With this concept, these waveguide problems become essentially 
two-dimensional, since for each mode the dependence on z (the direction of the wave- 
guide) is simply exponential. Furthermore, the number of modes is countably infinite. 
However, the situation is no longer so simple if it is possible for the electromagnetic 
wave to reach infinity in a direction perpendicular to the z-axis, as it can when the 
waveguide is open or when the waveguide has an imperfectly conducting wall. In the 
latter case, it is sometimes possible to use the idea of an impedance wall to save 
the situation. Otherwise, the modes lose at least some of their desirable properties, the 
excellent work of Marcuvitz and coworkers [1] in this connection notwithstanding. As 
an example of the open waveguide, consider the microstrip. Here the lowest mode, 
although exponentially decreasing in the direction of the z-axis, is unfortunately ex- 
ponentially increasing at infinity in any plane z = constant. This is readily under- 
standable and is intimately connected with the infinite length of the line assumed in 
the theoretical treatment [2]. On the other hand, for any physical waveguide or trans- 
mission line, the total length must be finite and consequently any meaningful solution 
of the physical problem must satisfy the usual radiation condition of Sommerfeld. For 
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this reason, the connection between the theoretical and the experimental aspects of 
the microstrip problem is by no means evident. Furthermore, it is apparently desirable 
to treat a problem of this type where the generator is not at infinity. No really simple 
problem of this type seems to have been treated in electromagnetic theory, and the only 
reasonable candidates seem to be the circular waveguide and the coaxial line. The 
advantage of the first is that it has one less dimension (the size of the inner conductor), 
while loosely the advantage of the second is that there is no mixing of the 7’ and the 
TM modes. Mathematically, the second problem is the simpler one of the two, and 
hence the imperfectly conducting coaxial line is treated in this paper. 

The precise geometry of the problem is shown in Fig. 1, where the imperfectly 
conducting outer conductor of the coaxial line is of finite thickness and runs from — © 
to + © along the z-axis while immediately inside a perfect conductor of zero thickness 
is fitted from — © to 0. To avoid meaningless complications, the inner conductor is 
assumed to be perfectly conducting from — © to + o. It is assumed that a TEM 
wave propagates from left to right; it is incident on the junction at x = 0. Thus there 
is effectively a generator at z = 0. This geometry is chosen because, in principle, at 
least, it possesses a closed solution in terms of quadratures. 

Since this problem is treated here only as an example of waveguide theory, there is 
no necessity of dealing with it in its full generality. Therefore, in the approximate theory 
of Secs. 5-10, it will be assumed that the free space wavelength is much larger than the 
transverse dimensions of the coaxial line in addition to the more specialized assumptions 
(5.1) and (5.2). 

This problem is also of considerable interest in connection with the so-called electro- 
magnetic shielding problem. However, in this connection, the assumptions (5.1) and 
(5.2) are too restrictive, and consequently, a more elaborate approximate theory has 
to be developed after a thorough understanding of the features of the more restrictive 
problem treated here has been acquired. 

Before studying the problem of the imperfectly conducting coaxial line, it is necessary 
to know some of the properties of the Green’s function to be used in the integral-equation 
formulation. This Green’s function is studied in the next section. 

2. Green’s function. The interest here is in the implications associated with the 
finite conductivity of the outer conductor. Since these implications do not depend 
critically on the choice of the dielectric constant of this imperfect conductor, this dielectric 
constant is to be identified with that of vacuum in the following consideration. If the 
conductor is also non-magnetic, then both the magnetic permeability and the dielectric 
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constant take those values appropriate for vacuum, say » and e. With the time dependence 
exp (— iw), Maxwell’s equations are 

V X E = wuH, 

V XH = (¢ — wO)E+ J, (2.1) 

V:-E=V-H=0O, 


Il 


where o is the conductivity of the imperfect outer conductor of the coaxial line and is 
replaced by 0 for air; the other symbols all have their usual meanings. If conventional 
cylindrical coordinates are used, all field quantities are independent of @ in this axially 


symmetric problem. Therefore 


With the continuity of 2, and H, as boundary conditions and the radiation condition at 
infinity, the electromagnetic field is completely determined by (2.1) provided that J 
is given. Let G,(r, ©) exp (itz) be the z-component of the electric field produced by the 
current source J = z 4(r — b) exp (7¢z), where z is the unit vector in the z-direction; 
it may then be verified that 


Go(b, §) = twul[k’E'Q, — (kh + ix’)n"Q.)", (2.3) 
where 
o, = HE GOHS eb) — HIGH Gb 9 
’ HS (a)HS? (&b) — HS? (ga) (Eb) ’ 
9, = | He?&e) Ho? (nb)Ho(ne) — Ho? (nb)Ho?(ne) | 
| neh (hk? + ix’) "HS? (0) 4" (nb)H5” (ne) — Ho” (nb)Ho” (ne) | (2.5) 
ye | He?) HS?’ (b)H2 (ne) — HE’ (nb)HS? (ne) 
nt E(k + ie) Hb) HE? (nd) HE (ne) — HE "(gd HS? (nd) | * 
g=(P-— 7)", n= K+ -— P)™, (2.6) 
and 
k? = wpe, «= woe. (2.7) 
Note that 
Go(b, $) = G(b, — 9), (2.8) 
and that 
Go(b, &) ~ —i(2kh? + tk’) we (2.9) 


as ¢— @, Because of (2.9), it is not meaningful to take the Fourier transform of G(b, ¢) 
with respect to ¢. Since G,(b, ¢) = O(k® — ¢°) asf — +k, it is possible to introduce 
an analog of the ordinary vector potential by defining 


g: (b, 2) = (2m)™ / dG, (b, 8) exp (itz), (2.10) 
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where 
G,(b, &) = (kb — &)G,(b, 9. (2.11) 
3. Formulation of the problem. The problem of the coaxial line may now be formu- 


lated in terms of a Wiener-Hopf integral equation. Let 7(z) be the z-component of the 
total surface current density on the infinitely thin perfectly conducting sleeve; the 


incident current density is exp (tkz). Thus as z — — @ 

i(z) ~ exp (ikz) + IT exp (—7ke), (3.1) 
where I is the coefficient of reflection with reference to the junction z = 0. With the 
g, of (2.10), 7(z) satisfies the integral equation 

d° 2 . pop , ra 

st+k dz’i(z’)q, (b, 2 — 2’) = 0 (3.2) 

dz Ja ; 


for z < 0. If & is assumed to have a small positive imaginary part which eventually 
approaches zero, then the Tourier transform of (3.2) in the form 


[ w| (2 — ix)ie || (2 + ik) (b,z — | = E,(b, z) (3.3) 
: Oz Oz 


—[(¢ — HO] [(e + HE, (6, O] = &.(b, O. (3.4) 


This equation is valid in a small strip | Im ¢ | < ¢, where G,(0, ¢) has no zero. The stand- 
ard Wiener-Hopf procedure calls for a factorization of G, in the form 


G(b, 9 = L.(H)/L_-(O), (3.5) 
where 
L.(¢) = exp § (2ri) ' [ - dt(t — 9) In G,(b, Op. (3.6) 
Note that 
L(g) = (L.(-—]°. (3.7) 
It follows from (3.4) and (3.5) that 
(kh? — PY(OLo = &.(b, OL_-(o). (3.8) 
This defines an integral function which is found to be a constant, say C. As a consequence 
of (3.1), the special case ¢ = k gives 
C = 2kL.(k). (3.9) 


From (3.8) the electric field is given by 


Eb, 2) = (2m) "C d¢{[L_(9)]~* exp (iéz). (3.10) 
In principle, the problem has been solved in closed form. 

4. The radiation pattern. The field far away from the origin is considered here. 
Since as £ > 0 


Ss 


Go(b, ) = —(twe)~'[b In (b/2)]€ [1 + OF Ind], (4.1) 
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H,(b, z) must be proportional to z~? as z > . This fact is of interest because it shows 
immediately that, as z— ©, an infinite number of “modes” must contribute significantly. 
The reason is that the combination of a finite number of modes must yield an exponential 
decay in the z-direction. This z~* decay may be understood as follows. Because the outer 
conductor of the coaxial line is not perfectly conducting, there is a coupling of the field 
inside the coaxial line with the radiation field outside of the coaxial line. Since the outside 
field decays as 2~” because of cylindrical symmetry, the field inside the coaxial line 
must also decay as 2~” 

Thus the term “radiation pattern’ is meaningful for the present problem. It follows 
from (3.8) and (3.9) that, as r, z > © with fixed @ = tan™'(r/z), 


H,(r,z) ~ A@(r? + 27)” exp [ik(r? + 2’)'”7], (4.2) 
where 
Siw L v (1) ° . ° (9)? \ (2 1)’ ‘ 
A(g) = Bie 42) 417.06 gin g) sin oH? (b®)H® (cb) — Hi? (b8)H®” (ca)] 
ac® L_(k cos $) (4.3) 
— k(k? + ix’) '@H§” ‘(ke sin ¢)|H5” (b®) Hs” (ce) — Hy” (b&)H," (c#)]}™", 
with 


& = (i? sin’ @ + ix’)'”. (4.4) 


5. Approximations for small skin depth. Although the exact solution has been 
obtained in principle, it is not very useful because of its complexity. The situation is 
greatly simplified if the coaxial line is only slightly leaky as when c — b is much larger 
than skin depth 6 = «~'v/2, or more precisely, if 


b> é (5.1) 
and 
exp [k(c — b) 2] > 1. (5.2) 
Under these assumptions, Q, of (2.5) is approximately 
Qs” = 4. (5.3) 
Furthermore when ke < 1, (5.3) leads to the following approximation of Go(b, ¢) 
Go?(b, §) = ino "[1 + Ck /8)")", (5.4) 
where 
C, = e'*“[kb In (b/a)]"’. (5.5) 


The physical meaning of (5.2) is as follows. When the coaxial line is lossless, the 
characteristic scale of length for the variation of the field vectors in the z-direction is 
k~*, which is very large. When the coaxial line is slightly leaky, the field inside the line 
tends to decrease when z increases, and this characteristic scale of length is the smaller 
one of k~' and xbk~" exp [x(e — b) 2]. If this characteristic length is much larger than 
the outer radius c, then the transverse problem is essentially separated from the longi- 
tudinal problem, except possibly near z = 0. Once the longitudinal and the transverse 
problems are separated, the Green’s function G, must admit simplification, and basically 
only the longitudinal problem need be considered. The condition (5.2) implies that the 
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frequency under consideration should not be too low in order that the thickness of the 


imperfect outer conductor be several times the skin depth. 
When i” is neglected compared with «*, (5.4) and (2.11) lead to 


Gb, ) = ie ike — 2) + C) — ey. (5.6) 
This may be factored by inspection to give 
LS(9) = (/0) ek + Ok + C,)'? + gy" (5.7a) 
and 
L? (8) = (/0) 7k — Ok + C,)'? — ¢). (5.7b) 
In particular, it follows from (3.9) and (3.8) respectively that 
CM se ag e771 — C,/4] (5.8) 
and 
8b, © = eo TL — C,/4fe re — ok + CY’? -— gy". (5.9) 


Note that &. is proportional to o~* as expected. 

The advantage of this approximation is that the factorization can be carried out 
explicitly. However, the exact Green’s function G,(b, ¢) as given by (2.3) has branch 
points at ¢ = + k, while the function G5’ (b, ¢) as given by (5.4) has no branch points 
there. Since it is precisely these branch points that determine the behavior of the field 
vectors for large z, this first approximation is not good enough. Since, as — — 0, the 


series expansion of the Hankel function yields that | H,’’(¢c)/H)5" (&) | = O(élné), 
a better approximation of Q, in (2.5) is 
Q;? = i(1 +P), (5.10) 
where 
P = 4H3) (ee [ne "(kb /x)? Ho” (&e) — Hy" (ée)]. (5.11) 


The term P is important only when £ is small. Therefore it is sufficient to use 


L? (kb) = LR, or C®? = Cc”. (5.12) 
Let 
M(¢) = Pll + C,(k/s)’)", (5.13) 
then (5.10) leads to 
G?(b, © = G!(b, O11 + MOT". (5.14) 
Let the quantity in the bracket of (5.14) be factored in the form 
1+ Mio) = 1 + N.O)/11 + N_(D], (5.15) 
where 
N.(¢) = —1 + exp 4 —(2mi)™’ , ; dt(t — ¢'In{il + M(o)} 


oe (5.16) 


= —(2Qri)™ | di(t — ¢)~'M(i). 


J —wowmie/2 











ag 
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It follows from (5.15) that 
L2(6) = LY (O[l + N.(9). (5.17) 


In the next few sections, the electromagnetic field at various points of the space 
will be studied. For simplicity all superscripts (1) and (2) will be omitted and only leading 
terms retained since, except at great distances from the junction at z = 0, the first 
approximation and the second approximation differ negligibly. 

6. Approximate field inside the coaxial line. In this section, the behavior of the 
field fora < r < bis considered. Fora < r < 3B, let 


cp — HP EDHEG) — HP GDH), “es 
= ~ HS? (€a)H§” (&b) — Ho”? (a) (Eb) ’ si 


so that 
&.(r, &) = &.(b, OSA, §). (6.2) 
Note that as — — 0, 
S,(r, ©) ~ In (r/a)/In (b/a). (6.3) 
Equation (6.3) implies that, when z > b, the transverse distribution of Z, is essentially 
In(r/a). 
The region a < r < b may be divided into various subregions I — V as follows 


Region I: — z >); 
gion III: z > b but kz | C, | « 1, and 
egion V: k | C,|z2> 1. 
Region IL is the region between regions III and I; and 
| Region IV is the region between regions III and V. 
The region II is the junction region and will not be studied further. 
| Let region I be considered first. In this region, F, is obviously very small. Thus, it is 
necessary to consider H,(r, z). It follows from (6.2) that the Fourier transform of H,(r, z 


’ 
Lc 
R 





1S 
Kar, £) = twet°*E,(b, $)(8/dr)S,(r, §). (6.4) 
If Res denotes “residue of’’, then in this region (3.1) holds with 
r = —[Res £°8,(b, §)]/[Res €°8.(6, )] = C,/4. (6.5) 
t=-k tak 
Thus 
H,(r,2) = Ts + 4C,e7**’). (6.6) 


For region IIT, it follows from (5.9), (6.2) and (6.3) that 


sid _, In(r/a) « se kz — 
Er, 2) = e**‘ko™ ” on exp | e sais . (6.7) 


In (6/a) xb In (ba) 


and 


H,(r,z) = ° e* exp je" a | (6.8) 


xb In (b a) 
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Therefore, the conventional theory of waveguides with a slightly lossy wall is valid 
provided that kz | C, | « 1 where C, is defined in (5.5). 
For the region V, the second approximation has to be used. The substitution of 
(5.17) into (3.8) gives 
F® (9) = F™ (ofl — N_(9). (6.9) 
When z is very large, only a small region in the neighborhood of ¢ = k can contribute 
to the inverse Fourier transform. Therefore, it is only necessary to consider the second 
term —F‘” (¢)N_(¢), since F“” (¢) is analytic at ¢ = k. With this notion, N_(¢) may be 


found explicitly near ¢ = k as follows: 

N_(g) = M(H), (6.10) 
since V,(¢) is analytic at ¢ = k. On the other hand, it follows from (5.13) that, near 
t= k, 

M(s) = S[k + C,) — &°)]4eck-*e"**”* In (€) exp [2ino(c — b)], (6.11) 
where 7) = x exp (iz/4). Consequently, for this region 
FP? (9) = —F (OM(0) (6.12) 


= 2k ao "(kl + C,/2) — &]7? ’c In (ec) exp [2inle — b)]. 
The substitution of (6.3) and (6.12) into (6.2) now yields 


In (r/a) 


’ c 27-3 —-1 19; 7 ’ ‘ 
4(7,2) =—-K«k xp [2 — b)| — a 
Er, 2) a Kh og exp | ino(e b)] thle”? (6.13) 
where 
F = | dge™[k(1 + C,/2) — 97h — °) In[(P — 2*)'e). (6.14) 


The integral in (6.14) may be evaluated by closing the contour of integration in the 
upper half plane. Thus 


F = —4miC{k*[1 — 2ni + 2 In (KC,c°’)] — QezC7k*[In (K’C,c°?) — in] 


2 (6.15) 
‘ait - | kC ane ; 
+ 4mik’e™* [ aze| Hs — | ot 


J0 ~ 


The first two terms on the right hand of (6.15) are never large enough to make a sig- 
nificant contribution so that they may be omitted. Thus 


- a al Calle ,, In(r/a) f°” re] kC er 
E,(r,z) = 4tcckk a fexp [2ino(e — b)}} ag | dite’ = =r) 2. (6.16) 


As z— © a stationary phase integration yields 
7 / 90.4 27.7-< kz as Ps . J 3 (pn ~ 
Efr, 2) ~ —32ix‘cb’k a7 'e'** {exp [2ino(e — b)]} In (b/a) In (r/a)z (6.17) 
and 


. 2 2 2 ike . / ] -—2 > , 
Hr, z) = 8ix “chk “e’™*{ exp [2ino(e — b)]} In (b/a) 72 =a (6.18) 
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Note that 2, behaves like z~*, and hence H, and FE, behave like z~*. It is interesting 
to note that the asymptotic phase velocity is that of vacuum. 

Finally, there is still region IV to be studied. Because of the form of Eq. (6.9), it is 
only necessary to add (6.7) and (6.16). Thus 


( 39/47, 
BE — - r/4 exp ( € kz ) 
4 + « ) : } 
xb In (b/a) (6.19) 


; pl ae re} AC ee i 
+ dieck='fexp [2in(e — DB} | dee" | i ‘| hate" nee. 


Jo 2 ) 
This is an approximate formula valid for all z so that z >> b. It gives explicitly the correc- 
tion to the simple theory (given by the first term) for waveguides with leaky walls. 

These explicit results varify in detail the discussion at the beginning of Sec. 4. 

7. Approximate field in the imperfect conductor. The region occupied by the im- 
perfect conductor may also be divided into five subregions, I‘ — V’, as given in the 
last section except that here b < r < ec. In region III’, the conventional waveguide 
theory with exponential dependence of z must hold again. Thus attention is concentrated 
on the region V’. 


The function that corresponds to the S, of (6.1) is now 
Sr, 2 = D(r)/D(b), (7.1) 
where 
1, (&) Hi (nr)H (ne) — HY? (nr) HS" ( a 
Dir) = F ‘ 7) (nc) 7 ") . (7.2) 
ne kk in’)? (ke) AG (a) HE (ne) — HY (nr) HS" (ne) 

Analogous to (6.2), it follows that, forb <r < ¢, 

&.(r, 2) = &.(b, OSAr, &). (7.3) 
For the present purpose, the appropriate approximate formula for S, is 
Sr, ©) = (b/r)'*fexp lin(r — b)|} (1 — exp [2in(e — r)] (7.4) 


2e~'**Lexp [2in(e — r)] — exp [2ino(e — b)] {xc(é/k)* In (ée)) 
when £// is small. With (7.3) written in the form 

&.(r, &) = 8:°(b, Of1 — M(O]SAr, O, (7.5) 
the electromagnetic field in region V’ is found to be 
Er, 2) = 4«“ebo'k*(b/r)'? In (b/ade"™” exp [in(r — b)](e'* ‘kf exp [Zine — r)] (7.6) 

— exp [2in(e — b)]}z~? — 8ixb In (b/a) exp [2in(e — b)]z*), 
and, with the second term in the parentheses neglected 
Hr, 2) = 4ix’chk-*(b/r)'? In (b/ade"™ {exp [ino(r — b)]}2- (7.7) 
-fexp [in(e — r)] + exp [in(e + r — 2b)]}. 

Here the second term in the brackets evidently represents a reflection at r = b. It is 
noted that (7.6) and (6.17) give the same value for E, at r = 6b, while (7.7) and (6.18) 
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give the same value for H, at r = b. This should be the case for a self-consistent scheme 


of approximation. 
A quanity of particular interest is the transverse distribution of the magnetic field 


in the outer conductor. For region IIT’, it is 


scoala ~1/8 —(r — b) —(2—-b-—r ‘“ 
| H,(r,z2) | ~r" | exp — + exp — Cea ba ny), (7.8) 
6 6 
while for region V’, it is 
elite oe .. call —(c — 1) —(e+r— 2b 
|H.(r,z2)|~r' | exp —s 5 dh exp =a tra 2b] (7.9) 


These distributions indicate clearly that, for regions III and III’, the electromagnetic 
energy effectively leaks from the coaxial line to the radiation field, while for regions V 
and V’, the electromagnetic energy effectively leaks from the radiation field into the 
coaxial line. This lends further support to the qualitative discussion of Sec. 4. For region 
IV’, the distribution of the magnetic field must be intermediate between that given 
by (7.8) and that of (7.9). These results are schematically sketched in Fig. 2. 

8. Approximate radiation field for large z. It remains to consider the situation 
outside of the coaxial line where the structure of the field is enormously complicated. 
In this section, attention is restricted to the extension of regions V and V’, namely the 
region V”: kz | C, | > 1andr >. This part of the calculation is therefore a continuation 
of that of the last section. In the next section the radiation pattern is to be found approxi- 
mately, as a continuation of Sec. 4. In Sec. 10, attention is turned to the extension 
III” of regions III and III’. There the relation with the lowest “mode” of the coaxial 
line is evident, and the present theory puts a limit on the range of validity of this “mode” 


solution. 
Forr > cand according to (2.3), the function that corresponds to S, and S, is given by 
4H,’ (ér) 
S¢,) = ————, (8.1 
s\") §7 rine D(b) tte 
and LE, is determined by 
E.(r, 2) = &.(b, OS3(r, $). (8.2) 


The situation here differs from the previous ones in that H,’’(ér), which is a factor in 























b c b c b c 
b << z << b/8 z~ rb/8 z>>Ab/8 


Fic, 2. Sketches of magnetic field distribution in the outer conductor, 
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(8.1), has branch points at ¢ = + kk. Furthermore, other contributions with branch 
points, from the determinant and from F(¢), have at least a factor &. For region V”, 
these other contributions may be neglected. Therefore, (5.9) is adequate for use in 
(8.2) with the result 


S.(r, &) = —2me' a ek bc!’ In (b/a) {exp [ino(e — b)]}H5" (er). (8.3) 


From the discontinuity of the Hankel function across the branch cut, the inverse Fourier 
transform gives 


EAr,2) = —2e'*a7'«*k *b?’c'? In (b/a) {exp [inole — nner *B(525), (8.4) 
where 
B(x) = [ s' dsJ,(s) exp (—is*x). (8.5) 


It does not seem possible to evaluate B(x) explicitly in terms of known functions. How- 


ever, it is easy to find that, as x — © 


iw te. (8.6) 
The substitution of (8.6) into (8.4) yields 
Er, 2) = 4e°°'0'*#°k 7 b?"c'”” In (b/a) { exp [ino(e — b)]}e"*27?. (8.7) 


This result is valid for z > b/ké and c < r « (z/k)'”*. Note that this field is independent 
of r, as must be the case for a far-sone radiation field. Also note that (8.7) is consistent 
with (7.6). 

It is not possible to find H,(r, z) directly from (8.7). To find H, , it is necessary to 
go back to (8.1), the reason being that 4dr nullifies the leading term in the present 
calculation. The determination of the magnetic field is not given here since it does not 
lead to anything interesting. 

9. The approximate radiation pattern. In this section, (4.3) is to be simplified 
under the assumptions (5.1) and (5.2). First, the functions L.(/) and L_ (k cos @) may 
be eliminated by (3.8), (3.9) and (5.9). The rest of the calculation is straightforward 
with the result 


A(o) = k’«*b'c' {exp [ingle — b]}[(l + C,)'? — cos ¢]' sing 
: a (9.1) 
ee a ke sin jee ss . 
{i —e' ‘ce sin’ ¢ In > O45. ‘xesin go | , 
where y 1.78107. Note that A(O) = 0, as it should. 
It is interesting to know the direction ¢, of the major lobe for this radiation pattern. 
This may be found as follows. From the last factor in (9.1), it is seen that 


, ke si ; 
xe sin’ dp» In Y on oe = O(1). (9.2) 
Therefore, in the vicinity of @) , A(@) is approximately proportional to 


*, +} “e 
A(¢) ~ const sin of — eke sin’ o _ é] ' (9.3) 
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The absolute value of this is 
yke sin 
2 
(9.4) 


+ «’c’ sin’ a(n = ny] 


This quantity on the right has a maximum at approximately 


| A(d) |? ~ const sin’ of — V2x«c sin’ ¢ In 


1 —1/2 : 
odo = (3 xe In “) (9.5) 
whence 
| —17.2,3/2,1/2 —(¢ — 0) | 1/271/2 ‘ 
A(g) | = «x kbc exp _— In (b/a)[{1 + 2°°°°)"". (9.6) 


10. The ‘‘mode.” Finally, consider the region III”’ defined by z > b, kz | C, | « 1 
and r > c. For this region (5.9) may be used for F(¢). Therefore, from (8.1) and (8.2) 
it follows that 

, 4H,’ (ér) 


S.(r, 0) = —e' no [kL + C,)'? — §] ; 
TUNC 
| Hs” (&e) H?? (ne) | 77 (10.1) 
| 
H;; (nb) » 2 
oF BOG) | 
; ) &¢ 0 if 
Ek + ik 7 
and hence 
Ar, §) = ek ER + Cy)’ — €)7"(b/0)'7 HH” (er) f exp [ingle — b)]} 
2 kc 2 k*s,\7' 
acl + In % _ tle . (10.2) 
T Z wKéC) 


In particular, when kr < 1, 3,(r, ¢) and hence H,(r, ¢) are proportional to 1/r. This 
fact agrees with the result given for the lowest ‘‘mode.’’ Secondly, when r increases, the 
phase of H;'’(r) cannot decrease near ¢ = k, and hence the phase of H, increases. This 
again agrees with the result for the lowest “‘mode’’ [2]. However, beyond that, the result 
for the lowest ‘‘mode’’ does not agree quantitatively with (10.2). In other words, for the 
present method of driving, the lowest “‘mode’’ does not dominate anywhere outside of 
the coaxial line. The ‘‘mode’’ picture therefore has very limited scope of application 
indeed. 

11. Summary and conclusions. The problem of the imperfectly conducting coaxial 
line driven as shown in Fig. 1 has been studied in some detail. The results for small skin 
depth may be summarized qualitatively as follows. The place where more detail may be 
found is written in square brackets. 

A. The reflection coefficient in the line at z = 0 is of the order 6/bln(b/a) [(6.5)]. 

B. For r < c and 6 K z € Ab/é the electromagnetic field decays exponentially 

[(6.7) and (6.8)]. 

C. In the dielectric inside the coaxial line, i.e., a < r < b, and for z > Ab/5, the 

longitudinal field decays as z~* but the transverse fields decay as 2” [(6.17) 


and (6.18)]. 











THE IMPERFECTLY CONDUCTING COAXIAL LINE 13 


In the outer conductor of the coaxial line and for z > \b/6, all field components 
decay as 2~* [(7.6) and (7.7)]. 
In the dielectric inside the coaxial line, the transverse distribution of each field 
component is essentially independent of z for all z > b[(6.19)]. 

. In the outer conductor the transverse distributions are not even approximately 
independent of z [(7.8), (7.9) and Fig. 2]. 

. Outside the coaxial line for z > \b/é and c < r « (z/k)'”, the field component 
E,, is approximately independent of r[(8.7)]. 

. The radiation pattern has a maximum in a direction slightly inclined from the 
direction of the coaxial line [(9.5) and (9.6)]. 
The ‘‘mode”’ picture has only qualitative meaning outside of the coaxial line 
[Sec. 10]. 

It is to be expected that many of these statements are true for any waveguide with 


a leaky wall and driven in a similar fashion. Let D, be the region enclosed by the wave- 
guide wall, D, be the region occupied by the waveguide wall, and D, = D, + D, . Also 
let L be a typical transverse dimension of the waveguide, then the following general- 


izations may be conjectured. 


The electromagnetic field decays exponentially for L « z<« AL /éin the region D; . 

Assume z > AL /6. Then for a TEM incident wave, the longitudinal field decays 

as z ° but the transverse fields decay as z-” in the region D, . For a non-TEM 

incident wave, all components decay as 2° 

or z > L/é and in the region D, , all field components decay as z~”. 

In the region D, , the transverse distribution of the field components are essen- 

tially the same for all z > L. 

‘’. The statement in £’ is not true for D, . 

. The radiation pattern has a maximum in a direction slightly inclined from the 
direction of the waveguide. (This seems to be true for both fast and slow waves.) 

’. The ‘‘mode”’ picture has very limited validity in general. 

The validity of these conjectures remains to be investigated. In particular, their 


correctness for a non-7EM incident wave should not be taken for granted. 


Acknowledgment. The author is indebted to Professors Ronald W. P. King and 
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BOOK REVIEWS 


The theory of storage. By P. A. P. Moran. Methuen Co., London, and John Wiley & 
Sons, Inc., New York, 1959. 111 pp. $2.50. 


This little book, the first in a new series of ‘‘Methuen Monographs on Applied Probability and 
Statistics’? presents on 105 pages a remarkably compact discussion of what is usually regarded as a 
rather rambling field: the mathematical theory of inventory control and water storage. In this country 
most attention has been focused on the first variant of this problem, the management of stocks of repair 
parts or other items with probabilistic demand so as to maximize some economic criterion of efficiency. 
This book offers a useful complement in that the dam interpretation of the problem, and those aspects 
of inventory control that bear most directly on water storage—on which the author is a leading authority 
—receive the most detailed treatment. 

An introductory chapter presents some probability concepts; but the reader should have some 
knowledge of the probability calculus beforehand. The chapter on the inventory problem emphasizes 
the calculation of probabilities of states, such as shortages, but pays little regard to the discovery and 
calculation of optimal policies. A useful and interesting feature is the systematic discussion of various 
inventory problems that do not involve stochastic processes but can be disposed of as one-period prob- 
lems. Segerdahl’s theory of insurance risk is also presented in outline. 

Throughout, the discussion is in terms of specific problems, and solutions are obtained for par- 
ticular distributions, typically those of the Gamma family. Continuous time problems are analyzed 
by going to the limit with discrete approximations. The Monte Carlo or simulation method is used to 
deal with the complicated problems that arise in the study of sequences of dams. 

In spite of its great variety of topics the book cannot achieve an adequate survey in its limited 
space. However, it is an excellent introduction on an elementary level to a field of operations research 
where much work is going on at the present. 





Martin BECKMANN 


Handbook of supersonic aerodynamics—mechanics of rarefied gases. By Samuel A. Schaaf 
and Lawrence Talbot. Johns Hopkins University Applied Physics Laboratory, 
Silver Spring, Maryland, 1959. ii + 85 pp. $1.25. 


The present section of the Handbook of Supersonic Aerodynamics is part of a series, the first one 
of which was issued in 1950. This section deals, as the title indicates, with rarefied gas flows. It opens 
with a general survey of rarefied flow regimes. Unfortunately, the original description of rarefied flow 
regimes as put forward by Tsien is still propagated. Following this is a useful chapter on free mole- 
cule flow calculations, in which the aerodynamic and heat transfer coefficients for some bodies of simple 
geometry are presented. The following chapter presents a discussion of slip flow which this reviewer 
considers to be somewhat dated in the light of more recent advances. This little compendium does have 
a considerable amount of information, much of it in graphical form, that may be useful to those in- 


terested in a cursory view of the field. 
XtoNALD F, PROBSTEIN 


German-English mathematics dictionary. Compiled and edited by Charles Hyman. 
Interlanguage Dictionaries Publishing Co., New York, 1960. 131 pp. $8.00. 


This useful dictionary contains more than 8500 entries from applied as well as pure mathematics. 
Extensive sampling did not reveal any incorrect entries. On the other hand, the temptation to increase 
the number of entries in a trivial manner has not always been resisted. For example, all of the entries 
“Axiom, Axiomatik, axiomatisch, axiomatisieren, Axiomatisierung, Axiomensystem’”’ could well have 
been omitted, and the space gained could have been used to list less trivial correspondences, for instance 
“Drall—moment of momentum”’ or ‘“‘Wirbelmoment—vortex strength.”’ 

W. PRAGER 


(Continued on p. 38) 
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ADDITION THEOREMS FOR SPHERICAL WAVE FUNCTIONS* 


BY 
SEYMOUR STEIN 
Sylvania Electronic Systems, Waltham, Mass. 


Abstract. Addition theorems are described for spherical vector wave functions, 
under both rotations and translations of the coordinate system. These functions are the 
characteristic solutions in spherical coordinates of the vector wave equation, such as 
occurs in electromagnetic problems. The vector wave function addition theorems are 
based on corresponding theorems for the spherical scalar wave functions. The latter 
are reviewed and discussed. 

1. Introduction. Boundary-value problems often involve several bodies or more 
than one important reference point. It can then be very convenient to be able to expand 
the field solutions in alternative forms, each form referring to a different specific co- 
ordinate set describing the same space. The connections between the alternate repre- 
sentations are provided by “‘addition theorems”, i.e., formulas for the expansion of the 
basis set of one representation in terms of the basis set of another. In general, 
the completeness of such basis sets indicates immediately the existence of these ex- 
pansions, although their specification may not be obvious. 

The results discussed here were motivated by a problem in electromagnetic scattering 
in which exactly such a need arose. The basis sets considered are the so-called spherical 
vector wave functions, i.e., characteristic vector field solutions of the vector wave equa- 
tion, based on a conventional separation of variables in spherical coordinates (R, 8, ¢). 
The set of vector wave functions can be described by using as potentials the set of 
scalar wave functions (the characteristic solutions to the scalar wave equation). Addition 
theorems for the latter are however available in the literature. In this paper, these are 
cited and briefly discussed, and extended to provide the desired general addition theorems 
for vector wave functions. The only other known calculations of these addition theorems, 
derived for specific low orders of vector wave functions by direct but laborious pro- 
cedures, and applied to some particular problems, appears in Ament [1]. 

2. Characterization of the vector wave functions. The electric and magnetic fields 
in a source-free homogeneous medium are divergenceless, and each satisfies a vector wave 
equation of form 

VxVxXxA-FA=0, (1) 
k being a constant for the medium. It is well-known (Stratton, [2]) that independent 
solutions of this equation can be constructed as 
M... = Vu. XR, 


17 x Man, (2) 


N.,.n — 
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with the added interrelation M,,, = 1/k V7 xX N,,, , where R is a position vector from 
the origin 0, and the potentials w,,,, are a corresponding complete set of solutions of the 
scalar wave equations 
Veut ku = 0. (3) 
In spherical coordinates, such a set of characteristic solutions are given by 
Umn = 2(kR)P™ (cos 6) exp (imd), —n < m<n, oA. < @. (4) 
Here z,(kR) is generically any spherical Bessel function, say as defined by Stratton 
3]. The associated Legendre function P” (cos @) will also be taken following Stratton 
[4], as 
, (1 a a)” 2 d' +m . A P 
Puig) = ao ssceie ig mi <n. (5) 
2'n! dx 
The latter differs by a (—)” factor from one common alternative definition. 
The specific forms for the vector wave functions are 
imP™ (cos @) 


; F . d 
exp (?m@) — i, exp (imd) — P” (cos 86) |, (6a) 
sin 60 ? diet ? dé 


M,.. = 2(6R)| i 


i 2, (kR) >m ; : ‘ : 
Nin = i LR n(n + 1)P% (cos 8) exp (2m@) 
ld x es d imP". (cos 6) ‘ 
—- (Rz,(kR)]| i, ex — p” Ss ; ho a ; (6 
LR dR [Rz,,( pI i exp (7m@) 716 P™ (cos 0) + i; a exp (ime) | ib) 


where i, , i, , i, are unit vectors in the directions of increasing R, 0, @ respectively. It 
is useful to note that forn = m = 0, we have M,, = Ny. = 0. 
We note that Maxwell's equations for harmonic time-dependence exp (— iw/), and 


for sourceless regions are 


V X E = lwuH , (7) 
é 
V XH = —weE, 
where w uve = kK’. It follows from Eq. (2) that to an £-field contribution of the form 


M,.,,, there is associated exactly an H-field term, (k/iwu)N,,,, ; similarly, to an £-field 
contribution N,,,, , there is associated an H-field term (k/iwu)M,,, . Furthermore, the 
detailed forms in Eqs. (6) show that M,,,, has no radial component, and hence all radial 
field components must be represented by the N,,,, solely. Thus, in fact, it is common to 
distinguish in electromagnetic theory two types of modes: the H-type in which only the 
magnetic field has a radial component, and the £-type in which only the electric field 
has a radial component. It is clear that H-type modes have E represented by the M,,, 
and H by the N,,,, , and the E-type modes vice-versa. For any general excitation, both 
types of modes may be present, and so the general field summations would be of the 
form 


B= > >Y [AnMa. + BasNun]; in 


n 1 m=—n 


H = (k/iwu) D> dS tAanNan + BanM nn]: (8b) 


Finally we note that if we make either a translation or rotation, or both, of the 
coordinates, a new set of scalar and vector wave functions can be analogously defined 
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with respect to the new coordinates. Since any one of the previous vector wave func- 
tions defines a perfectly valid vector field, it must be expandable over the new set. 
That is, under coordinates rotations and/or translations, there mus/ exist addition 


theorems of the form 
M,,, (2, 6, >) 
p oe C(u,v | m,n)M,,,(R’, 6’, ¢’) + Diu, v | m, n)N,,,(R’, 6’, ¢’)], (9a) 


’ 


N,,(R, 0, ¢) 
} > 7 C(u,v | m,n)N,,,(R’, 0’, ¢’) + Diu, v | m,n)M,,,(R’, 0’, o’)], (9b) 


where (/?’, 6’, ¢’) is the new set of spherical coordinates. In the above, we may note that 
the second equation follows directly from the first by utilizing the basic definitions in 
Eqs. (2), and noting additionally that the curl operator is invariant of the choice of 
coordinate system. We note also that the coefficients will of course contain implicitly 
the geometric parameters which relate the two coordinate systems. 

We will describe the representations for rotations and translations separately. A 
coordinate change involving both types will be a “‘multiplication”’ of the two operations. 

3. Rotation of coordinates. It is extremely simple to describe the addition theorems 
for M.,,, and N,,,, under coordinate rotations. 

We assume a second set of coordinates (R, 6’, ¢’) centered at 0, but with reference 
axes rotated with respect to the original (R, 6, ¢) system. Note that the radial coordinate 
is invariant under rotation. In Appendix 1, we characterize such a rotation, and state 
the known addition theorem for the spherical scalar wave functions, which is in the 


form* 
UnnA(R, 0,6) = Zz. B(u, m, n)ul,,(R, 6’, o’). (10) 


But the vector operator V is defined independently of the coordinate system, and due 
to the common origin, the same R is the position vector for both. It follows immediately 
by using Eq. (10) and the basic definitions that the vector wave functions are related by 


M,,,(2, 6,6) = Vum(R, 0,¢) XR = >> B(u, m, n)Mi(R, 6’, ¢) 


and this is the complete specification for coordinate rotations, of the general characteri- 
zation in Kq. (9). 

4. Translation of coordinate origin. Suppose a second origin of coordinates is taken 
at a point 0’, whose coordinates are (Ro , % , ¢0) with respect to the first. The set of 
spherical coordinates (R’, 6’, ¢’) is introduced with respect to 0’, such that the polar 
0, and the azimuth axis, ¢’ = 0, are respectively parallel to the corresponding 


axis, 6’ 
0 and @ = 0. This is then a rigid translation of the coordinate system. In 


axes, 6 
Appendix 2, we state the known addition theorem for the spherical scalar wave function 
under such translation, in the form 

*The prime notation on the umn is to emphasize the coordinate set with respect to which the po- 
tential is defined. We will utilize primes for this purpose throughout; under no circumstances should 


they be interpreted here as indicating differentiations. 
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u,(R, 0,¢) = p PR a(u,v | m, nul (R’, 0’, o’). (11) 


n=0 m=—n 

(We note from Appendix 2 that, depending on the value of F relative to R, , the u/,, 
may involve a spherical Bessel function 2{(kR) which is generically not of the same 
type as the z,(kR) involved in u,, . This distinction will not be important in the derivation 
below, which involves only the local properties of u/,,, , whichever Bessel function type 
is contained.) 

Since the grad operation is invariant of coordinate system, we can immediately 
write * 


M,,(R, 0,¢) = Vuw XR 
> D> alu,v| m, n)[Vur. X R]. (12) 


n=O m=-—n 
Our problem therefore is to expand terms of the form Vu/,, X R in terms of M/,, and 
N/,, vector wave functions. 
If R, is the position vector of the second origin, 0’, with respect to the first origin, 0, 


we can further write 


R=R,+R’. (13) 
But then 
Vuln XR = (Vuln X Ro) + (Vuln. X R’) (14) 
and we immediately recognize the last term as simply 
Vuln XR’ = Mi. (15) 


Hence it remains only to determine an expansion for Vu/,, XK Ro. 

The simplest approach to the latter has appeared to be to identify first the part 
which is represented by the N/,, functions, since this involves matching only the vector 
components in the ij direction. That is, using Eq. (6b), we seek an expansion of the form 


bead Dp 


-Vul. XRo= Dd Dd Aoil-N, , 


p=1 @=—p 
me Po Ag “or p(p + 1)P% (cos 6’) exp (igq¢’). (16) 
p=1 q=-D 
It may readily be ascertained that a fixed vector field of magnitude and direction 
equal to R, is represented at the general point (R’, 6’, ¢’) by components 
R,-i/ = R,[sin 0 sin 6’ cos (¢o — ¢’) + cos % cos 6’], 
R,-ii = R,[sin 0 cos 6’ cos @ — ¢’) — cos & sin 6’|, (17) 
R,-i; = R,fsin 4 sin @. — ¢’)]. 
Most of the derivation then consists in utilizing this result, along with trigonometric 
identities, and well-known Legendre and Bessel function recurrence relations. The 
algebraic manipulations are tedious but reproducible, and for the sake of brevity, the 
algebraic details will be omitted here. 


*In writing Eq. (12), we can recall explicitly that we are only interested in the form of expansion 
for vy > 1,| «| < »v; this is implicit throughout the remainder of the discussion. 
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For the radial part, comparison of the expansions on the two sides of Eq. (16) yields 


n= m= 0:if-Vusi XR, = 0, 


sia iia ._ i - 3 : 
n>1:%-Vu.xXR= ae iD i: \™ cos 6,N/,, 
_ = sin 6, exp (— ido) Nii41.2 (18) 


— = sin A Exp (ido)(n + m)(n — m + DN... 


This completely identifies the radial component, and hence all the possible N{,, which 
enter the expansion of Vui,, X Ro. 

If we next consider the if or if component of Vu/,, X Ro , subtracting off the corre- 
sponding component of the vector whose radial component is indicated on the RHS of 
Eq. (18), we can (after much manipulation) identify the remainder in terms of com- 
ponents of M/,, . The result may be summarized as follows: 


L=n=— ©: 
‘ , . kR, . — , , : 
V uso X Ro = kR, cos 65M4, + 9 sin 6, exp (—ido)M/, — kRysin 6 exp(ido)M’2,,, (19a) 


and forn > 1 
Vuln X Ro — — ee. Sim cos ON? — ~ sin 6) exp (— ido) Nii. 
_ : sin 0) exp (ido)(n + m)(n — m + DN...) 


= Mt + n+ Mae, n-m+I1y,, 
= M.. + 2n + - {208 af n Menn-1 + a +i mnt 


sin > exp (— ido) | -} :, _— 
+ 9 n Biaste~s + n + 1 Mas.042 
sin 65 exp (io) | + m)(n + m — 1) ,y,, 
ie < eeeel omcee ——-~weeemmnes 


a m—1,n—-1 
n 





(19b) 





(n—m+In-—m+%y, } 
n +1] m—i,n+1 ° 
We may note that Eq. (19a) is subsumed in Eq. 19b if it is assumed as a definition that 
M,,.-1 = 0 when n < 0 and that all terms in N,, are = 0 when n = 0; we will use these 
conventions in the last result below. 

We can refer back now to the general form of expansion in Eq. (9), and summarize 
Eqs. (11, 12, 14, 15 and 19). Further, making some obvious appropriate changes in 
certain of the summation indices, we find that the addition theorem under rigid co- 
ordinate translation has the form 


M,,(R, 6,¢) = D> Dd (Clu, »| m,n)Mi,.(R’, 6’, o’) 
ee (20a) 
+ D(u,v| m, n)Ni,(R’, 6’, ¢’)), 
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where 


kRy cos #.n + m+ I 


- (u,v | m,n 1) 
m+3 nti “ T 


C(u,v | m,n) = alu,v| m,n) + 


kR, cos 06.n — m 


(On — D a(u,v| m,n — 1) 


kR,sin 6, exp (—19) 
~ ? a(u,vy| m—1,n + 1) 


2(2n + 3)(n + 1) 
(20b) 
kR, sin 65 exp (—7i¢») 
+ : ? a(u,vy| m—1,n — 1) 
2(2n — 1)n 


4 kR, sin 0) exp (ido) (n + m + 2)(n + m + 1) 
7. ~ c 


: v(u,v | m l,n 1) 
2(2n + 3) n+ 1 ” + + 


kR, sin 6 exp (t6) (n — m — 1)(n — m) , 
— —— XP \’o we a(u,v| m+i,n —1), 
2(2n — 1) n 


ikR, cos 0 
D(yp,v | m,n) = - ma(u,v | m,n) 
“ n(n + 1) - 
ikR, sin 0) exp (—t@o) 
a a" —— c 


vu, v| m — 1,n) (20c) 
Qn(n + 1) “ 


tkR, sin Oy exp (to) ( 
2n(n + 1) 


In the above a(u, v | m,n) is taken to vanish whenever | m | > n (see, e.g., Eqs. (A2-4 and 


1+ m+ 1(n — malu,v! m+ 1,7). 


A2-2) of Appendix 2). 

It would of course be helpful if the expressions for C(u, » | m, n) and D(p, v | m, n) 
could be reduced, or could be shown to satisfy useful recurrence formulas. One might 
expect that use of formulas such as indicated in Appendix 2 might lead to such a reduc- 
tion. However, to date, although some suggestive results have been noted, no clear cut 
simplifications over the formulas above have been found. 

Acknowledgment. The assistance of Mr. George Foglesong, presently at MIT, in 
compiling some of the Appendix material, is gratefully acknowledged. 


APPENDIX I 
Addition Theorems for Scalar Spherical Wave Functions 
Under Coordinate Rotations 


These theorems are well known (especially in applications in quantum mechanics) 
as a special case of the study of the 3-dimensional rotation group. Various forms of the 
derivations, in largely non-group theoretical terms, appear in [5-9]. Some group theo- 
retical derivations appear in [10], along with general recurrence formulas for the co- 
efficients, and tables covering lower-order cases. We have found it convenient, however, 
to refer primarily to the forms common in quantum mechanics, such as presented by 
Edmonds [11]. In terms of the Euler angles of rotation, the addition theorem appears 


as 


Yin(9,¢) = >> Yim-(0’, ¢’) Do? (aBy), (Al-1) 
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with the relation to our notation [12] of 


2; Pa y]1v2 
Y,,.(0, ) = | Hi ') : i] P™ (cos 6) exp (im). (A1-2) 
The coordinates 6, ¢, are the angle coordinates of a point in space defined in the usual 
manner with respect to one set of cartesian axes, and 6’, ¢’ are similarly the coordinates 
of the same point with respect to a new set of axes, the latter set of axes being obtained 
by a rigid-body rotation of the first set of axes through the Euler angle a, 8, y. The con- 
vention used by Edmonds [13] is that the angles a, 8, y are positive in the sense of 
rotation of a right-handed screw in the right-handed frame of axes. The coefficients 
D\;,, (a@By) are the matrix elements of a unitary transformation, and are given by 


Dj, (aby) = [exp im’a] d,’,,(8)[exp imy], (Al-3) 


where 


(j YG — mt]? ; fe 
dg) = | G + m | (j — m’) " ( j J m )(? "7 
(j + m)!(j — m)! 7 \j—m’ —o@ o 


20+m'+m 2j;-2¢-—m’-—m 
j—-m'’-@ B . B 
(—) (cos 8) (sin 8) ; 


In terms of our form (Eq. 10), we note that both sides of Eq. (Al-1) involve surface 
harmonics of the same order, j, and hence it is a trivial step to introduce additional 
radial wave function factors z;(kR) on both sides to obtain a relation on total wave 
functions. Using our definition of the wave function (Eq. 4) we can rewrite Eq. (Al-1) 


(Al-4) 


in the form 
Un (Rk, 0,6) = 2. B(u, m, nun (R’, 0’, o’), (A1-5) 


where, using also Eq. (A1-2), 


min] (2 + wp)! (n — mI? 2 . 

B(u, m,n) = (—)” * menace tS ae D'S (aBy). Al-6 

“ ( (n — p)! (n + m)! Y ( ) 

Several useful sets of recurrence formulas are available for these coefficients, as special 
cases of general transformation laws [14-17]. but are omitted here, for brevity. 


APPENDIX II 
Addition Theorems for the Scalar Spherical Wave Functions 
Under Coordinate Translations 


1. A restatement of the Friedman-Russek result. A derivation of the translational 
addition theorem for scalar spherical wave functions has been given by Friedman and 
tussek [18]. However, they use a normalized Legendre function such that P[™ (cos 6) = 
P™ (cos @); and they appear to fail to note that their final statement of the addition 
theorem depends on the premise that | m + » | = | m| + |u| where m, yu are azimuthal 
indices of associated Legendre functions. Since the latter is not correct when m and u 
are of opposite sign, we cite here the basically correct form of their result, but give 
the correct values for the coefficients. In addition Friedman and Russek only identify 
these coefficients as the values of certain integrals over triple products of Legendre 
functions, the values of which in turn are found in the literature as rather complicated 
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sums. Instead, we point out that these coefficients are quite well known in the study of 
angular momentum in quantum mechanics, being in essence the vector-coupling co- 
efficients (or Wigner 3—j symbols) which describe the matrix elements of a unitary 
transformation between a product-state representation and a compound-state repre- 
sentation in a two-particle system. With this recognition, we are able to draw upon a 
large body of known information for useful computational forms and recurrence relations. 
Finally, we will correct below some errors of notation in the original results, as well as 
present a more useful form for one set of the results. 

To begin with, let us recall that our basic definition of the Legendre function (Eq. 5) 
can be regarded [19] as valid for negative m as well as positive, in the range —n <m< n 
(we have in fact so used P® in the text). 

With this in mind, the basic result derived by Friedman and Russek in [18] can be 
stated as follows: 

In an (R, 6, ¢) coordinate system, let the point 0’ at (Ry , 0 , do) be taken as the origin 
of a second coordinate system with coordinates (R’, 6’, ¢’), and oriented so that a rigid 
body translation [by the vector Ro(Ro , % , ¢o)| takes one system into the other. Then we 
have the addition theorem,* valid for y > 0, — v < » < »v, and where z, is any spherical 


cylinder function, 


z,(kR)P%(cos 6) exp (iu@) = 7” D> DS D> 4i"*7(2n + 1)(—)"a(u, m | p, », n) 


n=0 m=—-n p 


-z,(kr.)P5""(cos 6.) exp [i(u + m)¢,] (A2-1) 


f 
/-_ 
i 


‘jn(kr-)P, "(cos 6.) exp (— imé.)}) 


where 
r, = R' r-=R, 
10, = 0’ 6. = 04) when Rk’ > R,, 
d=¢' $= do 

and 


i. = @, @.= Or when & < £,. 
o> = do Ge = @ 
The sum over p is over all the values 


p=n+y, n+v— 2, n+v-— 4, 


(but no lower than | n — |), for which the coefficient a(u, m | p, v, n) is non-vanishing. 
The latter coefficient is defined by the expansion 
P™(cos @)P%(cos 6) = > 2 a(u m, | p, v, n)P?**(cos 8) (A2-2) 





*The introduction and use below of the notation 6, , @< , ds , @< serves to correct a rather obvious 
error in the Friedman-Russek paper, in the writing of their Eqs. (19) and (21); the error is obvious from 
their derivation which involves the parameters only in the quantities r cos y and ro cos yo. 
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which is known to exist in exactly such a form on the basis of spherical harmonic ex- 
pansion theorems, with n + vy > p > |n — v|. It is apparent from the orthonormality 
properties of the p% (cos 6) that a(u, m | p, v, n) can be determined from Eq. (A2-2) in 
terms of an integral over a special triple product of associated Legendre functions. Such 
integrals have been investigated in the literature of quantum mechanics, and Friedman 
and Russek cite one such evaluation, a rather complicated sum involving multitudinous 
factorials. 

Now, it is desirable to rewrite the result, Eq. (A2-1), in a more useful form, namely 
the form Eq. (11) in which we have used the addition theorem in the text 


u,(R, 0,¢) = yh Dd alu,v| m, n)ul,,(R’, 0’, ¢’). (A2-3) 
n=0 m=—n 
A comparison with Eq. (A2-1) shows that for R’ < Ry , we have exactly such a form 
where 


a(u,v| m,n) = 7 ’*"(2n + 1)(-)” ; Va(u, — m|p,v,n)us-"(Ro , 0 , 0). (A2-4) 


We note that in these equations for R’ < R, , the u/,, (R’, 6’, ¢’) contain explicitly a 
j.(kR’) dependence, no matter what the form of the spherical cylinder function in 
u,, (R, @,); and on the other hand, the uS~™ (Ro , % , do) contains exactly the same type 
of dependence z,(kR,) as is contained in u,, (R, 6, @) in the form z,(kR). 

However, for R’ > R, , the R’ variations in the form of Eq. (A2-1) go into the uf”, 
and in order to obtain a form like Eq. (A2-3), it is necessary to interchange the orders 
of summation. That such an interchanged form must be available is obvious, since the 
result is then exactly the one expected from the usual technique of expanding an arbitrary 
function of (R, 6, ¢) in terms of a new coordinate origin and new coordinate set. We 
accomplish the interchange by noting that in Eq. (A2-1), p eventually takes on all values 
between 0 and . Furthermore, the formal extension of the inner sum to 0 < p < @ 
can be made since the a(u, m | p, v, n) will vanish for all the added terms. We can also 
extend the sum on m to (— ©, @) since P® vanishes for all the added terms. But now, 
if we interchange the p- and n-summations, and substitute a new index ¢ = » + m, we 


obtain from Eq. (A2-1) 


ES feven + (-) als, tw lpm) 


<2 
n=0 t 


-z,(kr.)P5(cos 6.) exp (itd,) (A2-5) 


2,(kR)P*(cos 6) exp (tu) = 2°” _ 
-j,(kr-)P% ‘(cos 6.) exp [i(u — deal 


If we now interchange the letters p and n; write an m for ¢; and also note that P} = 0 
whenever | ¢ | > p, we can rewrite the above as 
2s af . 
2,(kR)P%(cos 0) exp (ive) = 7" DE DY D4 2"*"(2p + 1)(—)" “alu, m — w| 0, p, ») 
n=O m=—-n p=0 \ 


-2,(kr,)P"(cos 6.) exp (im@,) (A2-6) 


-j,(kr-)P%- "(cos 6.) exp [t(u — méap 
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It is readily confirmed that the coefficient in this last equation also vanishes un- 
less, at least, n + »v > p > | n — v|, so that it is even more in a form similar to Eq. 
(A2-1). By now comparing with Eq. (A2-3), we see that for R’ > R, , the result corre- 


sponding to Eq. (A2-4) is 
a(u,v!| m,n) = 27 °*"(-)”"* > (2p + l)a(u, m — w!n,v, pur "(Ry , 00 , 0), (A2-7) 
= 


where now u,-” contains explicitly the factor j,(kR,) while the more general dependence 
on position has a z,(kR’) of exactly the same cylinder function type as z,(kR). 

Although this result appears to be different from the form in Eq. (A2-4), one expects 
on the basis of continuity of the two expansions across the surface R’ = R, that they 
should be equivalent. That this is indeed so, and hence that Eq. (A2-4) truly represents 
the required coefficient for all R’, is most easily shown by referring to the more general 
representations of a(u, m | p, v, n) coefficients as “vector-coupling coefficients’ or 
“Wigner 3—j symbols’. The equivalence can then be shown to follow from symmetry 
properties of the Wigner 3—j symbols [20]. 

Again special recurrence relations are available [21], while others can be derived 
from Eq. (A2-2) by using well-known Legendre function recurrence formulas. 
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SPIN MATRIX EXPONENTIALS AND TRANSMISSION MATRICES* 


BY 
L. YOUNG 


Stanford Research Institute 


Abstract. The three Pauli spin matrices o,(¢ = 1, 2, 3) occur in the mechanical, 
especially quantum mechanical, theory of rotation in three-dimensional space. The three 
spin matrix exponentials are here defined as exp (¢;2), where x is the independent vari- 
able. Transmission matrices can be expressed in terms of spin matrix exponentials, 
thereby permitting a more systematic treatment of transmission line circuits. 

Introduction. In the design of quarter-wave transformers, it has hitherto always 
been assumed that the guide wavelength is independent of position along the line. 
This is so, for instance, for TEM modes; or for TE», modes in rectangular waveguide 
where the wide or ‘a’ dimension is kept constant. Such transformers, having guide 
wavelength independent of position, are called homogeneous transformers [1]. The 
first exact design formulas for ideal homogeneous quarter-wave transformers were 
given by Collin [2], who considered up to four sections. (The junction of two trans- 
mission lines when junction discontinuities are neglected, is called an “‘ideal transformer’. 
This is analogous to two perfectly coupled coils of turns ratio (Z,/Z,)'’* and having 
infinite inductance.) The first complete synthesis procedure was given by Riblet [3]. 
The author later computed extensive numerical tables [4], which have been checked 
out experimentally on numerous occasions. 

tiblet’s synthesis procedure [3] is based on Richards’ transformation [5] and Richards’ 
theorem [6], and thereby depends on the commensurability of all transmission line 
sections in the circuit. The homogeneous quarter-wave transformer has also been used 
as a prototype circuit in the design of direct-coupled-cavity filters [7]. 

It has been shown that the performance of single-section quarter-wave transformers 
can always be improved by going from a homogeneous to an inhomogeneous design [8]. 
The analysis of inhomogeneous transformers of more than one section has only recently 
been undertaken [9], and the purpose of this paper is to present the mathematical tools 
which were developed for this purpose. A separate paper will deal with the design 
considerations and numerical results for multi-section inhomogeneous quarter-wave 
transformers [21]. 

Spin matrix exponentials. With line-lengths no longer commensurable, a more 
general formulation than is possible by Richards’ transformation is required. For a 
systematic and compact treatment of transmission matrices, we shall employ the three 
Pauli spin matrices [10], which may be represented as follows: 


*Received May 9, 1960. This paper was written while the author was at Westinghouse Electric 
Corporation; it is based on part of a dissertation for the degree of Doctor of Engineering at the Johns 


Hopkins University, Baltimore, Md., 1959. 
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0 ] 
Gg, = 
1 0 
fo —1 
a aed j I (1) 
] 0 
( 
1 0| 
G3 = | 
0 -l) 
They anti-commute among themselves, and their squares are equal to unity: 
G10, = —6q0; = jos 
0203 = —0302 = Jo; (2) 
ey “ese = jo2) 
o=-0,=0;=I], (3) 


where J is the unit matrix (idemfactor). We define the ‘“‘spin matrix exponentials” by 


Ex) = exp (ao;) = I cosh x + o; sinhx (¢ = 1, 2,3) (4) 
and their derivatives by 
. De 
E'(x) = © EX(2) = ¢; exp (xo;) 
dz 


(5) 
= ¢;E;(x) = I sinh x + o; coshz (¢ = 1,2, 3)) 
They do not commute among themselves unless 7 is the same, and then they commute 


and behave like ordinary exponentials: 


E<xE<y) = Ex + y). (6) 
Also, 

E(2)o; = o;E;(x) (7) 
but 

E(a2)o; = 0;E;(—2), t ¥ j. (8) 
Similarly, 

Et(@E(y) = E(x t+ y), (9) 

EX(z)o; = o,Ei(a), (10) 

EX(z)o; = —o,Ei(—2), 1}. (11) 


Spin matrix exponentials occur in the mechanical (especially quantum mechanical) 
theory of rotation in three-dimensional space [11, 12]. The matrix product 


Q = £(#)p,(#),(#), (12) 


where 6, ¢, ¥ are the three Eulerian angles, yields a matrix whose four elements are the 
Cayley-Klein parameters [11]. At this point, we anticipate in order to complete the 
analogy: it will be seen later that Eq. (12) is like the transfer matrix of a single-section 
ideal transformer. This is not altogether surprising, since geometrical analogies have 
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been developed before [13, 14, 15], and projective charts have been used for the numerical 
solution of transmission line problems. The spinor theory of two-ports and its geometrical 
interpretation has also been discussed in general terms by Payne [16]. 

The Pauli spin matrices, together with the unit matrix, can be used to express any 
matrix [17] in the form of a quaternion, but the resulting more general form does not 
have the simplicity of the spin matrix exponentials which are well-suited for the analytic 
description and solution of inhomogeneous transformers. 

The transfer matrix. There is no uniform terminology for the transformation 
matrices [18-20] which are used to analyze two-ports. The matrix defined by 


“) = 7||, (13) 
b, ibe 


where a, , }, are the incident and reflected wave amplitudes at the input, and a, , b» 
those at the output (Fig. 1), will be referred to as the transfer matrix. As only lossless 
two-ports are here considered, wave amplitudes may be defined in terms of power by 


a |° = power flow in the forward direction 


(i.e. towards the load) g (14) 
b |? = power flow in the backward direction! 
(i.e. towards the generator) J 


The transmission coefficient, 7’, between two reference planes is the same (in both 
phase and amplitude) when going from left to right as when going from right to left. 
The transfer matrix can then be written [18] 





I _r: 
tf." mf (15) 
QM pp _ Pils 
gy d Yi 


where 7’ is the (unique) transmission coefficient, IT, is the reflection coefficient seen 
at the input (on the left) when a matched load is placed at the output (on the right), 
and I, is similarly defined for the output side. 
We may, by appropriate choice of reference planes, let 
r,= -r,=TPF (Gay). (16) 


From energy considerations: 
1. Let b, = 0 in Fig. 1. Then 


IrP+ {iT l? =1. (17) 
ao) -+ i | —— 
b, TWO -PORT Sinn b> 
INPUT OUTPUT 


(2) «2 (2) 


Fia. 1. Defining the transfer matrix. 
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2. Let a. = Oin Fig. 1. Then 
(7 +r!] =1. (18) 


Comparing Eqs. (17) and (i8), we infer that, under the condition (16), the trans- 
mission coefficient vector is parallel to the reflection coefficient vector [22]: 


T parallelto TI. (19) 


Now further choose the reference planes so that 


lr = real. (20) 
The transfer matrix then reduces to: 
25 Se x 
T=5| e (21) 
rF 4 


The diagonal and anti-diagonal parts of T. From Eq. (15) it follows that for a 
reflectionless two-port, T must be a diagonal matrix. Let 





pa jiu Py2| (22) 
[T, T2) 
Two useful concepts are 
iT, 0 
Di(T) = |" | (23) 
0 T'22) 
= diagonal part of T 
eo Ba 
Ag(T) = | ma (24) 
Fox Of 
= anti-diagonal part of T. 
For zero reflection, 
Ag(T) = 0. (25) 


The ideal transformer. The magnitude of the reflection coefficient at an ideal 
transformer (Fig. 2) is 


r| = ie — 5 (26) 
Z,+ Z, ”" 


Select both reference planes to be coincident in the plane of the junction itself, and 
choose T to be real. This satisfies Eqs. (16) and (20) and therefore leads to the form (21) 
for the transfer matrix. T is now determined except for sign. In agreement with Ref. [4], 
we (arbitrarily) pick the positive sign and let 


Ze or. Z (27 


: = oe 
Z2 + 2; 
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POSITION OF REFERENCE 
PLANES (REAL [) INDEPENDENT 
OF FREQUENCY. 


a 


J 











| 

| 

! 
—eoro 
ond 


INPUT OUTPUT 
SIDE Zz, Ze SIDE 
Z, ,Z2: CHARACTERISTIC IMPEDANCES. 
Fic. 2. Ideal transformer. 
Defining the “junction VSWR” by 
— * ‘ 
V = Z, (28) 
and the “log ratio’ by 
a=tinV (29) 
the transfer matrix (21) of the ideal transformer reduces to 
T = E,(a) (30) 
which is the first spin matrix exponential as a function of the log ratio, a, of the junction. 
Also T = tanh a, (31) 
and since T is real, 
T = (1 — I’)'” = secha. (32) 


Length of transmission line. A section of transmission line of electrical length @ 
radians (Fig. 3) has a transfer matrix 


T = £,(98) (33) 


which is the third spin matrix exponential of imaginary argument j@. 
The ABCD matrix. The ABCD matrix of an ideal transformer is 


A = E;(@) (34) 
and that of a length of transmission line is 
A = E,(j@). (35) 


The transition between transfer and ABCD matrices has thus been effected merely 
by interchanging suffixes of the first and third spin matrix exponentials. 








Fic. 3. Section of transmission line. 
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Conclusion. The spin matrix exponentials exp (o; x), where z is either j0(@ = elec- 
trical line length) or the log ratio of a transformer, represent the transformation matrices 
of line sections and transformer steps. They are useful in the treatment of transmission 
line transformers, particularly inhomogeneous transformers, where the line lengths are 
incommensurable and Richards’ transformation then does not apply [21]. 
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ON COMPLETE SOLUTIONS FOR FRICTIONLESS 
EXTRUSION IN PLANE STRAIN* 
BY 
J. M. ALEXANDER 
Department of Mechanical Engineering, Imperial College 


Consideration is given to the extension of partial slip line field upper bound solutions 
to give the true yield point load when the material is constrained. In particular, the 
problem of frictionless extrusion is studied and it is shown that it is possible to extend 
only one of three available partial solutions. If it is not possible to extend the available 
partial solution, the use of discontinuous stress fields leads to lower bound solutions, 
and examples of this technique are given. , 

1. Introduction. It is well known that all slip line field solutions constitute upper 
bound solutions, being partial incomplete solutions as described by Bishop (1953). 
This is because the slip line field is always associated with a kinematically admissible 
velocity field, although it is not generally associated with a statically admissible stress 
field, as discussed by Prager (1959, p. 118), for the extrusion problem. Bishop (1953) 
has proposed techniques for extending such partial stress solutions into neighbouring 
rigid regions, with a view to demonstrating whether or not a statically admissible stress 
field exists for the partial solution considered. If such a statically admissible stress field 
does exist, then the partial solution provides a lower bound also, giving in fact the true 
yield point load for the problem concerned. It is a consideration of the possibility of 
extending partial solutions to give true yield point loads for the problem of frictionless 
extrusion which forms the subject of this paper. 

2. Available partial solutions, 3:1 ratio. For the 3: 1thickness ratio of sheet extrusion 
through a square-edged die there exist three possible solutions, as shown in Fig. 1. 
\ll three solutions were proposed by Hill (1948) and in discussing the solutions shown 
in Fig. 1b and Ic he pointed out that the existence of dead metal regions would require 
the satisfying of certain conditions by the friction called into play on the surface of 
the dead metal. 

Now the solution shown in Fig. 1(a) has the lowest yield point load, and since all 
three partial solutions are differing upper bound solutions, it should be possible to 
extend only that given by Fig. 1(a). This may not in fact be possible since there may 
be a further partial solution giving a still lower upper bound. 

3. Complete solution, 3:1 ratio. A possible statically admissible stress field is shown 
in Fig. 2, for the partial solution of Fig. 1(a). Also shown in this figure is the stress plane 
showing the cycloidal traces of the pole of the Mohr’s circle, using Prager’s (1953) 
geometrical representation. The broken line 11 — 6° is a principal stress trajectory, 
and the stresses are transmitted across this boundary and supported on an infinite 
number of rectangular elements, via triangular elements of the type shown inset, the 
faces of the rectangular elements being normal and parallel to the extrusion axis. This 
concept was used by Bishop in the paper already referred to. The principal stress o, 


*Received May 13, 1960; revised manuscript received June 27, 1960 
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on any triangular element can be supported by two direct stresses p, and p, both equal 
to a, , So that no shear stresses are invoked. The rigid material to the left of the principal 
stress trajectory will not yield provided that | p, — p. | < 2k. The maximum value 
of this stress difference occurs for the shaded element and is of amount (¢,),; — (¢,):0° 
where (o,);,; and (¢,),9- are the principal stresses acting normal to the principal stress 
trajectory at the points 11 and 10’ respectively. The broken line 0A is a stress dis- 
continuity, as suggested by Bishop, the stress states on each side of it being as indicated 
in the diagram. Since the infinity of columns transmit no shear stress, the boundary, 
equilibrium, and yield conditions are everywhere satisfied and the solution is a statically 
admissible stress field constituting a lower bound. Since it is also an upper bound, the 
actual yield point load will be associated with it. 

Considering the solution shown in Fig. 1(b), it can be shown immediately that an 
associated statically admissible stress field does not exist, as follows. The angle 8 is 
less than 2/4, and since the shear stress on the vertical die face must be made zero to 
comply with the requirements of the stress boundary conditions, the rigid dead metal 
vertex at 0 cannot support the boundary shear stresses postulated without yielding. 
That this is so can be seen from the analysis by Hill (1954) (case iv), in which it is proved 
that the angle 8 cannot be less than 7/4. 

Considering the solution shown in Fig. 1(¢), an attempt to extend the partial solution 
is shown in Fig. 3. On the left hand side of the principal stress trajectory 10’ — 6’, 
the maximum stress difference (in the rectangular element A) is less than 2k. On the 
right hand side of the envelope of slip lines 9 — 10’, an equilibrium stress state on each 
triangular element can be as shown in the inset diagram, from which it is seen that 
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DP, = p —k, p. = p + k, where p is the mean stress at any point on the envelope. Thus 
at the point 9, p, = po — k, and at the point 10’, po = pio» + k, as shown in the stress 
plane. Since the Mohr’s circles for these two points are certainly not coincident, the 
element B is greatly overstressed in that p, — p2 > 2k. This is true for the whole of 
the region 10’ — R — Q — 9 of the rigid dead metal, so that this region would yield 
under the distribution of stress applied to it. It is interesting to note that Hill’s (1954) 
analysis is not applicable to the whole of the vertex 9 — 10’ — S, since the problem 
of the variation of stress on the face of such a vertex is not there considered. This ex- 
tension of Bishop’s infinity of rectangular elements permits the consideration of the 
possible yielding of such assumed rigid domains. 

4. Extrusion ratios greater than 3:1. In Fig. 4 is shown the extended partial solution 
giving the true yield point load for extrusion ratios greater than 3:1, together with 
the stress plane, using the methods just described. The complete range of validity of 
this solution has not been investigated. 

5. Other extrusion ratios. In Fig. 5 is shown a statically admissible stress field for 
an extrusion ratio of 2:1. This solution is interesting in that it has a dead metal zone, 
even with a frictionless interface between container and material. There is a field of 
constant stress in the dead metal region, the stress state being as shown in the inset 
element. It has not been found possible as yet to find complete solutions for other ratios. 

6. Discussion. It has been demonstrated that the techniques suggested by Bishop 
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for extending partial solutions can be applied to problems of contained deformation. 
In this connection, it is the author’s opinion that Bishop’s extension of Hill’s field sug- 
gested in his paper (Bishop, 1953) is unacceptable because it violates the stress boundary 
conditions at the frictionless container—material interface. Problems of contained defor- 
mation with friction at such interfaces require special consideration, particularly con- 
cerning whether or not relative motion exists at such boundaries. The formulation of 
lower bound theorems for such cases constitutes a difficult problem. 

It may often happen that none of the available partial solutions can be extended to 
give a statically admissible stress field, and hence the true solution. Under these cir- 
cumstances a field of stress discontinuities in statical equilibrium may be found which 
would give a lower bound solution. 

As examples, two discontinuous stress fields for the 3:1 ratio frictionless extrusion 
problem are shown in Figs. 6 and 7, in which the broken lines are the stress discon- 
tinuities separating regions of constant stress. The magnitudes of the constant stresses 
in each region are shown by the inset elements whose sides are parallel with the principal 
planes in the particular region concerned. In the stress plane of each diagram are shown 
the Mohr’s circles for the states of stress in each of the regions, together with their 
‘poles’, as defined by Prager (1953). The slip lines of both these fields intersect all bound- 
aries and the axis of symmetry at angle 7/4, so that they are statically admissible, 

lower bound solutions. The field shown in Fig. 6 leads to a mean extrusion 
pressure of 8/3k, whilst that of Fig. 7 gives 10/3k. (The field shown in Fig. 7 was derived 
from that suggested by Shield and Drucker (1953) for the indentation problem). 

For this particular problem it has been possible to show that the upper bound slip 
line field of Fig. 2 gives, in fact, the true vield point load of the deformation, namely 
a mean extrusion pressure of 24/3 (2 + 7) = 10.284/3k. Had this not been known, 
however, taking the mean extrusion pressure midway between the best available upper 
and lower bounds would have resulted in an error of about 13% only. 


and therefore 
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Upper bound solutions for many problems of metal working have been proposed 
by Johnson (1959); the lower bound techniques proposed here form a complementary 
method of attacking such problems. The possibility of extending the upper bound 
methods to deal with three dimensional problems has been suggested previously, Alex- 
ander (1959); the use of discontinuous stress fields has been shown by Shield and Drucker 
(1953) to be suitable for application to three dimensional cases for finding lower bounds. 
Thus methods are available for bounding yield point loads for three-dimensional metal- 
working problems, which is a field requiring further study. 

Acknowledgment. The author is indebted to Professor Hugh Ford, Professor of 
Applied Mechanics, Imperial College of Science and Technology, for many helpful and 
stimulating discussions relating to this work. 
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The numerical treatment of differential equations. By L. Collatz. Third Edition. Springer- 
Verlag, Berlin, Gottingen, Heidelberg, 1960. xv + 568 pp. $23.56. 


Though labelled ‘third edition’, this is essentially an English translation of the second German 
edition (see this QUARTERLY, vol. 13, p. 348). The references have been brought up to date, many 
improvements have been made, and the number of examples has been increased. 
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Operations research and systems engineering. Edited by C. D. Flagle, W. H. Huggins, 
and R. H. Ray. The Johns Hopkins Press, Baltimore, 1960. x + 889 pp. $14.50. 


Part I (Perspectives) is concerned with the purpose and the historical development of operations 
research and systems engineering. Part II (Methodologies) is primarily devoted to the mathematical 
techniques. It contains chapters on simplified models (E. Naddor), basic statistics (A. J. Duncan), 
statistical quality control (A. J. Duncan), digital computers (W. C. Gore), inventory systems (E. 
Naddor), linear programming (V. V. McRae), queuing theory (C. D. Flagle), simulation techniques 
(C. D. Flagle), theory of games (E. Naddor), symbolic logic (W. E. Cushen), design of experiments 
(W. G. Cochran), human engineering (A. Chapanis), information theory (W. C. Gore), flow-graph 
representation of systems (W. H. Huggins), system dynamics (W. H. Huggins), and feedback and 
stability (N. H. Chosky). In Part III, the use of these principles and methods in operations research 
and systems engineering is illustrated by case studies. 
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Progress in solid mechanics. Volume I. Edited by I. N. Sneddon and R. Hill. North- 
Holland Publishing Co., Amsterdam, 1960. xii + 448 pp. $15.50. 


This is the first volume of a new series, which is to be primarily devoted to the ‘‘basic principles and 
mathematical techniques of continuum mechanics, in all its aspects, together with experimental work 
of a fundamental kind’’, As the space available for this review is totally inadequate for a critical discussion 
of the eight articles, the following highly condensed table of contents will have to indicate the character 
of the volume. 
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ON THE THEORY OF THE PLASTIC POTENTIAL* 
By HANS ZIEGLER, Ziirich 


1. Introduction. The theory of the plastic potential, proposed by v. Mises [1], 
connects the yield condition and the flow rule of a plastic solid. In its generalized form 
established by Prager [2], the theory may be stated as follows: 

Let q, (k = 1, 2, --- ,n) be the generalized strains, i.e. a set of coordinates specifying 
the deformation of the plastic body in the sense of analytical mechanics, and let q; , 
g; denote respectively the elastic and the plastic components of q, . If the work done on 
an infinitesimal increment of strain is given by 


dA = Q, dq. ’ (1.1) 


where the summation convention has been applied, the Q, are the generalized forces 
in the sense of analytical mechanics or, according to Prager, the generalized stresses 
corresponding to the coordinates q, . 

The generalized states of strain and stress may be depicted by points with position 
vectors q, Q in Euclidean n-space R,, , or by the vectors q, Q themselves. Defining the 
scalar product of two vectors in R, by 


Q-q= AG ; (1.2) 
we obtain for the work (1.1) the representation 
dA = Q-dq. (1.3) 


Consider the solid in an arbitrary stage of the deformation process, and let Q denote 
the actual state of stress. Any state of stress Q* which can be reached from Q without 
plastic flow will be called nonplastic. All of the Q* constitute the nonplastic domain 
R,, in R, , and any stress increment of the type Q* — Q will be denoted as nonplastic’). 
The yield limit is a hypersurface in R,, , defined by those nonplastic states of stress the 
infinitesimal changes of which are not exclusively nonplastic. 

The theory of the plastic potential firstly stipulates that the elastic strains q° follow 
from the stresses Q according to the laws of elasticity. The plastic strain increment 
dq’ corresponding to given values of the stress Q and the stress increment dQ is zero 
for any state of stress lying in the nonplastic domain but not at the yield limit. For 
states of stress at the yield limit, dq” may be different from zero; it is secondly stipulated 
that 


(Q* — Q)-dq’ < 0 (1.4) 
for any nonplastic stress increment Q* — Q, and that 
dQ-dq’ > 0. (1.5) 


In virtue of (1.4) the nonplastic domain RF, is convex. 


*Received June 17, 1960. 
'In Prager’s or Koiter’s terminology [3], 2, is the elastic domain and Q* an allowable state of stress. 
The notations used here emphasize the fact that ideally plastic as well as hardening solids are considered. 
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2. Formulation of problem. So far, the theory of the plastic potential is but a 
hypothesis. This has been emphasized already by v. Mises; besides, his version was 
restricted to a volume element with strains e;; and stresses o;; . The generalization of 
Section 1 is due to Prager; the special form (1.4), (1.5) does not refer to a potential 
function and was given by Drucker [4]. 

Several authors have tried to provide a basis for the theory (see [3], p. 180). Bishop 
and Hill [5] have shown that (1.4) follows for the element of a polycrystalline aggregate 
from plausible assumptions concerning the behaviour of the single crystals. Drucker [4] 
has based the theory on a postulate regarding the work done on a prestressed element 
by an external agency in a cycle of application and removal. The author [6] has suggested 
a generalization of Onsager’s theory of irreversible processes [7] to nonlinear cases, 
thus providing a thermodynamic basis for the theory of the plastic potential. 

Either one of these approaches is based in its turn on certain postulates some of 
which are open to criticism [8]. Thus, it seemed worthwhile for once to reject any kind 
of postulate and to limit the scope of the investigation to a purely mathematica! proof 
for the theory in Prager’s generalized form, based on the sole assumption that it holds 
in v. Mises’ sense for an element of volume. This proof has been provided by the author, 
in [8] for rigid-plastic materials and in [9] for elastic-plastic solids. The present paper 
gives a condensed version of the two articles and at the same time provides a simpli- 
fication of the proof. 

3. The rigid-plastic element. In a rigid-plastic solid, q’ = 0 and hence q = q’. 
The local states of strain and stress are given by e¢;; , ¢;; respectively and may be depicted 
by the vectors e, s in Euclidean 9-space R, . Defining the scalar product in R, by 


s-e = 0 5 ;€;; 9 (3.1) 


we obtain for the work per unit volume done on an infinitesimal increment of strain 
(3.2) 


dA = o;; de;; = 8-de. 
In virtue of the symmetry of the strain and stress tensors, the vectors e and s actually 
lie in a linear subspace R, of R, . Here, the nonplastic domain may be defined according 
to the rules of Section 1. If we assume that the theory of the plastic potential is valid 
in v. Mises’ sense, it follows from (1.4) and (1.5) that 


(s* — s):de = (o* — a;;) de;; < O (3.3) 
for any nonplastic stress increment s* — s, and that 
ds-de = do;; de;; > 0. (3.4) 


These relations also hold (with the equality sign) for states of stress below the yield 
limit. In virtue of (3.3) the nonplastic domain is convex. 

4. The rigid-plastic solid. The states of strain and stress of the whole body B are 
given by the functions e;;(2,), o;;(2,) respectively, either one of them depending on 
the coordinates x,. These functions may be represented by vectors E, S in function 
space F. Let the scalar product in F be defined by the volume integral 


7 


S-E = | o4;€;; AV (4.1) 


/B 


extended over the entire body B. This definition is admissible (see [10]), since it satisfies 


the commutative law, S-E = E-S, the distributive law, S-(E + E’) = S-E+ S-E’, the 
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rule S-O = 0 for multiplication by zero and the associative law for scalars, (aS)-E = a 
(S-E). Besides, the definition provides F with a positive definite metric. On account of 
(4.1), the work done in an infinitesimal increment of strain is 


ye i) 6 = OR. (4.2) 
B 


The representation considered here is not restricted to kinematically admissible 
states of strain nor to statically admissible states of stress; it holds for any states for 
which the integrals in (4.1), (4.2) exist. The nonplastic domain F in F may again be 
defined according to the rules of Section 1. In general, certain elements of B reach their 
local yield limit for states of stress S still inside the yield limit of the whole body. Plastic 
flow sets in when a sufficiently large domain of B has become plastic. At this stage, the 
state of stress in any element of B either lies below or on its local yield limit. Since none 
of the elements undergoes plastic flow in a nonplastic stress increment of the entire 


body, it follows from (3.3) that 


(S* — S)-dE = [ Gt~ alt, dv <6 (4.3) 
7B 


for any nonplastic stress increment S* — S. Also, on account of (3.4), 


. 


dS-dE = | da;; de;,; dV > 0. (4.4) 
7B 


Hence, the theory of the plastic potential, if valid for the element, likewise applies to 
the body as a whole. Incidentally, relations (4.3) and (4.4) also hold (with equality 
sign) for states of stress below the yield limit. In virtue of (4.3), also the nonplastic 
domain F is convex. 

5. Generalized strains and stresses. Section 4 provides a basis for the rigorous 
treatment of a rigid-plastic body. In numerous cases, however, one is compelled to 
simplify the problem by introducing generalized strains and stresses as defined in Section 
1. This means necessarily that only states of strain are considered which can be de- 
scribed by a finite set of parameters q, (k = 1, 2, --- , n). Such a reduction of the degree 
of freedom can be realized by introducing additional constraints; in certain cases also 
the elimination of originally existing constraints may result in a simplification (see 
examples in [8] and [9]). It is obvious that this process is only justified as long as it 
does not appreciably modify the actual state of strain. 

The limitation considered here implies that only states of strain E are considered 
which belong to a certain subspace F,, of F. It does not involve a similar restriction with 
respect to the states of stress. Let E“ denote the state of strain in function space F 
corresponding to the generalized strain gq, = 1, q; = 0 (¢ ¥ k). The subspace F,, then is 
defined by the states of strain 


E = ¢E"’. (5.1) 


° k . ° . 
It is reasonable to assume that the vectors E“ are linearly independent. In this event, 
there is a one-to-one correspondence between the vectors E in F,, and q in #,, . On account 
of (4.2), (5.1) and (1.1), the work done by the stress S on an infinitesimal strain increment 


dE belonging to F,, is 
dA = S-dE = S-E” dm = Q. dy . (5.2) 
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Hence, the generalized stresses are given by the scalar products 
Q, = S-E™. (5.3) 
It follows immediately that any state of stress S in F is represented in FR, by a unique 
vector Q which, conversely, is the image of an infinity of vectors S. 
The straight section connecting the points in F with position vectors S’,S”’ is given by 
S = S’ + a(S” — 9’), O<a<l). (5.4) 


If Q, Q’, Q” are the vectors in R, corresponding to S, S’, S’’, we obtain from (5.3) and 
(5.4) 

QO. = Qi + a(Q’ — QD, <4 «6s 0). (5.5) 
Thus, the image in R, of the section (5.4) is the straight section connecting the points 
with position vectors Q’, Q”. 

Since the nonplastic domain F is convex, any vector S* — S representing a nonplastic 
stress increment lies entirely in F. Its image in R, is the vector Q* — Q which, according 
to the definition of the nonplastic domain, lies entirely in R, . Conversely, any vector 
Q* — Q in R, is the image of at least one vector S* — S connecting two points of F. 
Since F is convex, S* — §S lies entirely in F; hence, any vector in R, can be considered 
the image of a nonplastic stress increment. 

From (5.2) and (4.3) we obtain 


(Q* — Q)-dq = (S* — S)-dE < 0, (5.6) 
where Q* — Q is an arbitrary nonplastic stress increment. Likewise, on account of 
(5.2) and (4.4), 

dQ-dq = dS-dE > 0. (5.7) 


This is the proof that the theory of the plastic potential, if valid for the element, like- 
wise applies to the treatment of the entire body in generalized strains and stresses. In 
virtue of (5.6), the nonplastic domain R, is convex. 

6. The elastic-plastic solid. In Section 1, no rules have been specified for the de- 
composition of the strain into its elastic and plastic components. In the case of a volume 
element, however, the decomposition is straightforward, provided the strains are 
sufficiently small’. Let us postulate that in the expression 

dA = s-de* + s-de’ (6.1) 
following from (3.2) the first product represents the infinitesimal increase of elastic 
strain energy per unit volume, while the second one represents the work dissipated in 
the infinitesimal strain increment de. Then, e° is the strain corresponding to the stress 
s according to the law of elasticity, and e” is the permanent strain still present after 


removal of the stress. 
In order to define a similar decomposition for the finite body, let us postulate that 


also in 


dA = S-dE‘ + S-dE’ (6.2) 





2For a few critical observations, see [9]. 
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the products are respectively the increment of elastic strain energy and the work of 
dissipation. Thus, E° is built up from the local elastic strains e* and E” from the local 
plastic strains e”. Here, E* may be different from the state of strain corresponding to the 
given loads if the body were elastic; also, E” is not necessarily the state of strain after 
removal of the loads. On the other hand, the definition (6.2) implies that the work of 


dissipation is 
dA’ = s-de” = o;; dé; (6.3) 


per unit volume and 
dA? = SdE = [ 04; dé, aV (6.4) 
JB 


for the entire body. 

Let us now replace (3.2) by (6.3) and (4.2) by (6.4). Assuming that relations (3.3) 
and (3.4) hold for the plastic strain increment de” instead of de, and retracing the demon- 
stration of Section 4, we arrive at (4.3) and (4.4) with (dé; instead of de;; and) dE” 
instead of dE. 

For the representation in generalized coordinates, the situation is similar. Let us 
postulate that the products in 


dA = Q-dq° + Q-dq’ (6.5) 


are respectively equal to the increment of elastic strain energy and the work of dissi- 
pation. Then (5.2) can be replaced by 

dA” = S-dE’ = Q-dq’, (6.6) 
and the remainder of Section 5 leads to (5.6) and (5.7) with (dE” instead of dE and) 
dq” instead of dq. 

Thus, the results proved in Sections 4 and 5 also hold for the plastic strains in elastic- 
plastic solids, provided the elastic and plastic strain components are defined by means 
of strain energy and dissipation work’. 

7. Observations. In an elastic-plastic solid, the elements which reach their local 
vield limit are not surrounded by rigid material. Hence, plastic flow sets in as soon as the 
first elements become plastic. It follows that under otherwise identical circumstances 
the nonplastic domain of the elastic-plastic solid is usually distinct from the one of the 
rigid-plastic body. Moreover, any plastic flow is apt to modify the nonplastic domain. 
Hence, the yield limit of the elastic-plastic solid undergoes a continuous transformation 
as the plastic flow proceeds towards collapse. On the other hand, the yield limit of a 
non-hardening rigid-plastic solid is always the same. 

In order to determine the nonplastic domain R, in an arbitrary stage of the defor- 
mation process, it is necessary to consider the actual state of stress S. From the general- 
ized stresses Q alone no information concerning the yield limit is available. In the case 
of a non-hardening rigid-plastic body, the yield limit must be determined only once. 
Once it is known, the theory of the plastic potential may be applied in the sense of 
Section 1, and the problem can be treated henceforth in generalized strains and stresses. 
In fact, this is the reason why Prager’s version of the theory offers an essential simpli- 


2As, e.g., in [6]. 
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fication. In the case of an elastic-plastic solid, however, one is compelled to keep track 
of the continuous changes in shape of the yield limit. Since this requires that the actual 
state of stress S be pursued throughout the deformation process, the use of generalized 
strains and stresses does not seem to offer any advantages here, except, of course, in 


limit analysis. 
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MOMENTS OF THE GENERALIZED RAYLEIGH DISTRIBUTION* 


BY 
J. H. PARK, JR. 
General Mills, Inc., Mechanical Division 
Minneapolis, Minnesota 


I. Introduction. Gaussian processes are of considerable interest in problems involv- 
ing random noise. Also of interest is the Rayleigh distribution which arises in work on 
radar, the detection of signals in noise, ete. [1, 2]. The generalized Rayleigh process 
promises to be of interest in the future especially when signals in noise are thought of 
to exist in a finite dimensional Hilbert space [3, 4, 5]. The generalized Rayleigh process 
was defined and some of its properties were investigated by K. 8S. Miller, et al. [6]. The 
purpose of this paper is to investigate the moments of the generalized Rayleigh dis- 
tribution. 

Let X, , X, , --+ Xy be independent Gaussian random variables with means Z, , 
Zy respectively and equal variances of one. A generalized Rayleigh random 
(also referred to as a non-central chi-square variable) is defined as 


y?= > Xx? 1.1) 


ro 


variable , Y, 


and the density function of Y, denoted g(y), is given by [Ref. 6, Eq. 1.6] 
[yoy/yo)*”? exp [(yo + y*)/2Mww-2/2(yoy) for y > 0 (1.2) 


gly) = 
0 for y £0, 


where 
N 
y= D8 (1.3) 
i=1 
and J,(x) is the modified Bessel function of the first kind. g(y) is called the generalized 
Rayleigh distribution. In this paper expressions for the moments about zero of g(y) 
and several interesting properties of these moments will be derived. 
It does not complicate the problem to consider non-integer moments. Therefore, 
the ath moment. of g(y) is given by 


Ss) 


M.(N, w) = [ yo) dy, (1.4) 


where ‘‘a’’ is any real number. (However, as will be seen later, the above integral exists 
only fora > — N hence a can be any real number greater than — N.) Whenever only 
integral moments are considered the subscript » will be used. 


The important results are: 


*Received April 27, 1960; revised manuscript received September 19, 1960. This work was supported 
by the U.S. Air Foree under Contract No. AF 33(616)-3374 at the Radiation Laboratory, The Johns 


Hopkins University. 








46 J. H. PARK, JR. [Vol. XTX, No. 1 


) The power series expression for 1/7,(N, yo), 


ne rT (N jz : 
MAN, ys) = 2°7 exp (—v2/2) TEE OE ay, (1.5) 





2) The “closed form” expression for J/, (N, ny 


M.AN, yo) = 2°” exp (—y5/2) ee ae M((N + a)/2, N/2, yo/2], (1.6) 


where M (with no subscript) is the confluent hypergeometric function (in the notation 


of Jahnke and Ende, [7, p. 275)). 
) The asymptotic expressions for 1/, (N, y), 


T om | a — 2) _ 
Yo) ~ vf + a(N ta— 2) + - = 2N h = Se + "as | 


M,N, 1? 


2y 2! (2yé)? 
as Yo ©, (1.7) 
and 
M.(N, yo) ~ N°’ exp (ays/2N) as N- @. (1.8) 


(4) The recursion formulas 





M.sAN, yo) = (N + a@)MN, yo) + ui] MV, Yo) + Yo aM .(N, Yo) | (1.9) 


dyo 
and 
M.-AN, yo) = 4 exp (—y5/2)(y5/2)°*~?” | si exp (x) **-°"M[(22)""*] dx 
Jo 
for a>2—N. (1.10) 
(5) The upper bounds on negative integer moments 

M.. & | 1/u% |" n=1,2,---N-—2 (1.11) 
and 

M., = |1/WV -—7n) |” n=1,2,---N. (1.12) 


II. General expressions for the moments of the generalized Rayleigh distribution. 
We first obtain a simple expression for the moments as defined by Eq. (1.4). Substituting 
(1.2) in (1.4) and abbreviating M, (N, yo.) by M, we obtain 


(2 2 /a\ a 2a+N)/2 j > V6 / . 
M,=y exp (— Yo/2) | y exp (—y°/2)T cw-2)/2(yyo) dy. (2.1) 
v¥WVJ 
The integral in the above expression diverges if a S — N. Therefore, whenever a 
is used it will denote a real number greater than — N and n will denote an integer greater 


than — N. Substituting in (2.1) the equivalent power series for J.y+2)/2 (yyo) [8, p. 163], 


we obtain 


oy N+2r+a-1 2 sé 
' , ’ Te exp (— 7 eA : 
_ |) , eT lL - : $ 9) 9” 
Mi =2 exp (—yo/2) | > TIN + 27/2 (yo/2)*" dy. (2.2 


[(J aS 


Interchanging summation and integration in the above, we note that the integral is a 


gamma function, hence Eq. (1.5). 
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The sum in (1.5) ean, recalling that Br(B) = T(B + 1), be rewritten, 
> Ir+Wet +_a)/2] (2, 2)" _ Ii(W+4)/2) 
r! Tir + (N/2)] Yo/# T'(N/2) 
{1+ 54 Tov ++ aan + 2» asay}. (2.3) 


The term in brackets in (2.3) is the power series for the confluent hypergeometric func- 
tion, M[(N + a)/2, N/2, y3/2], and will sometimes be abbreviated by M. Therefore, 
using (2.3) in (1.5) we obtain the closed form expression (1.6). 





III. Asymptotic expressions. In this section asymptotic expressions for M,(N, y,) 
as a function of y with a and N fixed and as a function of N with a and y fixed will be 
derived. 

First consider the case when a and WN are fixed. /,(N, yo) is given by (1.6) in terms 
of the hypergeometric function and therefore, using the asymptotic expression for Mf 
given on p. 275 of [7], we obtain (1.7), the asymptotic expression for 17,(N, yo) as yo >. 

To obtain the asymptotic expression for M, as N — © we use the power series 
expression for J/, given in (1.5) where the series has been rewritten as shown in (2.3). 
The product in the right hand side of (2.3) can be written 


TT yar rr ar eee > (N+2k+a)N 
I] (N + a+ 2k)/(N + 2k) = [1 + (@/N)] IT w + aN +h 


viii 2ka 
-u+@my fit - aa 3ta |: 6-0 


=0 








It is clear from (3.1) that this product is asymptotic to [1 + (a/N)]". Hence from (2.3) 


and (1.5) we obtain 





P(N + 0@)/2) 0 paar ; 
MAN, Yo) ~ ~~ Tv72)_ 2°” exp (ayi/2N) as No. (3.2) 
But 
TUN + @)/2) par joya/ ar 
ry N/2)"? as N—@ (3.3) 


and hence Eq. (1.8). 


IV. Recursion formulas. Recursion formulas can be easily derived for moments of 
order a + 2 in terms of the moment of order a and its derivative with respect to yp . 
Since M,(N, yo) = 1, moments of even integer order are easy to compute from this 
formula. A recursion formula for moments of order a — 2 is also obtained in terms of 
an integral involving the moment of order a. 

As a preliminary step, consider Eq. (1.5) where x has been used in place of y?/2. 


After differentiating with respect to x and rearranging terms we obtain 


dM. a ~ I(r + (N + a)/2] 
2M, + St) = 2-20 rr + (V2) 





; (4.1) 








But 





Mu. = 22%? exp (—2) >s p+ W ro ei + (N+ 4)/2) (4.2) 
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Therefore 


Mis2 =(N+a@M,+ 2o( M, + aM). (4.3) 
4 1 


If we put (4.3) in terms of y) we obtain the recursion formula given in (1.9). In particular 
f £ ' 


since M, = | 


M_=N+ ym, 
M, = N(N + 2) + 2y(N + 2) + yi, (4.4) 
ete. 


To obtain the recursion formulas for decreasing moments we multiply M,(2) by 


exp (x)x‘****’? and integrate term by term from 0 to #?/2 resulting in 
auo?/2 
| exp (z)z°"*°-?°M,[(22)"""] dz 
‘ '4a-2)/2qa/2 Ta Tir — 1 (N €1)/2], 2 . 
iia Zz. — = bs : " (y/2)’. (4.5) 
, 0 r! Tir + (N/2)] , 
But 
' rir—1+(N a) /2 oe 
M,-2 = 2°~?” exp (—y3/2) D> | a + 3/3 (yo/2) (4.6) 


r! T'[r + (N/2)] 


which together with (4.5) gives Eq. (1.10). It can be shown that the integration in 
(4.5) exists only ifa > 2 — N. 
When a = 0 (1.10) becomes 
M_AN, yo) = 3 exp (—y;/2)(y5/2)°-*” | exp (aja'*~*?? de (4.7) 
and for a few particular values of N we have 
M_.(4, yo) = 1/y, 
M_.(6, yo) = {1 — (2 Yo)| Yo 
and 
M_(8, yo) = [1 — (4/ys) + (8/y5)]/ys- (4.8) 


V. Upper bounds on negative integer moments. ‘Two very simple expressions which 
are upper bounds on negative integer moments can easily be derived. Let n be a negative 


even integer, say 7 — 2m, m = 1, 2, --- then 1.5 can be written 
| . er " 
9) , am — 7° /9) 2 /a\" 9 T a, .. Oh & 
M _om 2 exp (— Yo/2) p i a (Yo/ 2) I] 2/(N + 27 2k). (5.1) 
7 u ° h 1 
Assuming (V/2) — m > 0 the right side of (5.1) is not decreased by letting r be zero 


in each term of the product resulting in 
Mom S&S 1/(N — 2m)”. (5.2) 


Alternately if we assume (NV/2) — 1 2 m we can replace m for (N/2) — 1 in the product 
in Eq. (5.1) and not decrease the right hand side. Therefore we have 
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M om < 2°” exp (—y5/2) > = 3 (y,/2)". (5.3) 
on ow a or +m)! ~ 
The sum in the right hand side of (5.3) can be written 
I —- r! ' a“) Wee 
(yo /2) = (2/y,)” (yo/2)’. 5. 
> r! y (r + m)! Yo) »» r! Yo) <) (5.4) 


By filling in some positive terms in the sum on the right hand side of (5.4) it becomes 


an exponential and we obtain 


M2. 3 1/y"- (5.5) 
By Schwartz’s inequality 
Miu & M2Mx. (5.6) 
In particular if m = 0 and k is a negative integer, i.e., k = — n,n = 1, 2, --- then 
Mr. & Mu. (5.7) 


Combining (5.7) with (5.5) and (5.2) we obtain (1.11) and (1.12) respectively. 


BIBLIOGRAPHY 


1. S. O. Rice, Mathematical analysis of random noise, BSTJ (1) 24, 46-156 (1945) 
2. J. L. Lawson and G. E. Uhlenbeck, Threshold signals, MeGraw-Hill Book Co., Inc., 1950 
3. W.H. Huggins, Representation and analysis of signals, Part 1, The Johns Hopkins University Report 


No. AFCRC Tr. 57-357 (1957) 

1. J. H. Park, Jr. and E. M. Glaser, The extraction of waveform information by a delay line filter technique, 
1957 IRE Wesson Convention Record, (2) 1, 171-184 

5. David Middleton, A note on the estimation of signal waveform IRE Trans. on Information Theory, 
IT-5, 86-89 (1959 

6. K. 8. Miller, R. I. Bernstein, and L. E. Blumenson, Generalized Rayleigh processes, Quart. Appl. 

Math. 16, 137-145 (1958) 

Eugene Jahnke and Fritz Ende, Tables of functions, Dover Publications, New York 

8. N. W. McLachlan, Bessel functions for engineers, Oxford University Press, London, 1934 











50 BOOK REVIEWS [Vol. XIX, No. I 
BOOK REVIEWS 


(Continued from p. 38) 


Continuous distributions of dislocations, by B. A. Bilby (Burgers vector and torsion tensor—Shape, 
lattice and dislocation deformations—Zero lattice pure strain—Dislocation density and stress—General- 
ized spaces). 


Asymmetric problems of the theory of elasticity for a semi-infinite solid and a thick plate, by R. Muki 
(Solution of the equations of equilibrium by Hankel transforms—Solution of the thermo-elastic equations 
by Hankel transforms—Stresses in a semi-infinite elastic solid under the compressive action of a rigid 
body—Stresses in a semi-infinite elastic solid with a tangential load on its surface—Thermal stresses 
in a semi-infinite elastic solid and a thick plate under steady distribution of temperature). 


W. PRAGER 


Dynamic programming and Markov processes. By Ronald A. Howard. Technology Press 


of M. I. T., and John Wiley & Sons, New York and London, 1960. viii + 136 pp. $5.75. 


Consider a physical system S represented at any time ¢ by a state vector x(t). The classical description 
of the unfolding of the system overtime uses an equation of the form z(t) = F(z(s), s < t), where F is a 
prescribed operation upon the function z(s) for s < ¢. In certain simple cases, this reduces to the usual 
vector differential equation dz/dt = g(x), x(0) = c. 

For a variety of reasons, it is sometimes preferable to renounce a deterministic description and to 
introduce stochastic variables. If we take x(t) to be a vector whose i-th component is now the probability 
that the system is in state 7 at time ¢, and allow only discrete values of time, we can in many cases 
describe the behavior of the system over time quite simply by means of the equation z(t + 1) = Az(t). 
Here A = (a;;), 7,7 = 1, 2, +++ , N, is a transition matrix whose element a;; is the probability that a 
system in state 7 at time / will be found in state 7 at time ¢ + 1. Processes of this type are called Markov 
processes and are fundamental in modern mathematical physics. 

So far we have assumed that the observer plays no role in the process. Let us now assume that in 
some fashion or other the observer has the power to choose the transition matrix A at each stage of the 
process. We call a process of this type a Markovian decision process. It is a special, and quite important, 
type of dynamic programming processes; cf. Chapter XI of Dynamic Programming, Princeton Univer- 
sity Press, 1957. 

Let us suppose that at any stage of the process, we have a choice of one of a set of matrices, A(q) 
(a;:;(q)). Associated with each choice of g and initial state 7 is an expected single-stage return );(q). 
We wish to determine a sequence of choices which will maximize the expected return from n stages of the 
process. Denoting the maximum expected return from an n-stage process by f;(n), the principle of 
optimality yields the functional equation 


N 
f,(n) = max | bg) + > a.,(@fi(n — 1) 
q i=1 

In this form, the determination of optimal policies and the maximum returns is easily accomplished 
by means of digital computers; cf. S. Dreyfus, J. Oper. Soc. of Great Britian, 1958. Problems leading 
to similar equations, resolved in similar fashion, arise in the study of equipment replacement and in 
continuous form in the “optimal inventory’’ problem; see Chapter Five of the book mentioned above 
and K. D. Arrow, 8. Karlin, and H. Scarf, Studies in the Mathematical Theory of Inventory and Production, 

Stanford University Press, 1959. 
As in the case of the ordinary Markov process, a question of great significance is that of determining 
the asymptotic behavior asn — ~. It is reasonable to suspect, from the nature of the underlying decision 
process, that a certain steady-state behavior exists as n — ©. This can be established in a number 


of cases. 


(Continued on p. 62) 


HIGH FREQUENCY VIBRATIONS OF CRYSTAL PLATES* 
BY 
R. D. MINDLIN 


Columbia University 


1. Introduction. In this paper, Cauchy’s [1] two-dimensional equations of coupled 
flexural and extensional motion of crystal plates are extended to the next higher order 
of approximation so as to accommodate the frequencies of the two lowest thickness- 
shear modes. The equations obtained are also extensions of previous equations [2] in 
which flexure and the same thickness-shear modes were included, but coupling with 
extensional modes was omitted. The new equations are deduced from the three-dimen- 
sional equations of linear elasticity by a procedure based on the series expansion methods 
of Cauchy [1] and Poisson [3] and the variational method of Kirchhoff [4]. Comparison 
of the appropriate solution of the resulting equations with Ekstein’s [5] solution of the 
three-dimensional equations, for an infinite plate, reveals close agreement between the 
two frequency spectra, over the extended range of frequencies, for all five branches of 
the spectrum of the two-dimensional equations. This indicates that solutions of these 
equations, for bounded plates, will give reliable results, over the extended range of 
frequencies, since the modes of bounded plates, in that range, are composed essentially 
of coupled overtones of the first. five modes of vibration of an infinite plate. 

In addition to the derivation of the approximate equations, theorems of uniqueness 
and orthogonality are established and some general conclusions are drawn regarding 
solutions in rectangular coordinates and vibrations of rectangular plates. 

2. Expansion in power series. The plate is referred to rectangular coordinates 
zx; (¢ = 1, 2, 3) with x, and x; in the middle plane and the faces at x, = + h. The com- 
ponents of displacement u; (j = 1, 2, 3) are expanded in power series of the thickness- 


coordinate: 


" (n) 0) (1) 2. (2) 
u; = >> xu" =u, +2; +2; °°, (1) 
n=0 
where the u;” are functions of zx, , x; and the time, ¢t, only. 
Stress-equations of motion. The series expression for u; is substituted in 


[ (Ti3.3 bo pu; es) bu; dV =0 (2) 
y 


which is obtained from the variational equation of motion, deduced from Hamilton’s 
principle [6]. In (2), the integration is over the volume, V, of the plate; the 7;; are the 
components of stress; p is the density; and the summation convention for repeated 
indices is employed, as is the comma notation for differentiation with respect to the 
z, and the time tf. 

When the integration with respect to x. , from — h to h, is performed in (2), the 


result is 


*Received September 26, 1960. 
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[ > (7: pale? +P? — 9X H gas) bu\” dA = 0, (3) 
A n=0 


2 m=O 


where A is the area of the plate and 
ah 
rry(n) nryy y(n) ney jh 
T” = | 27 ,, dt,,  F” = [x2T,,]*, 
h 


h 


2h"*"*"/(m +n + 1), m+n even (4) 
Hon. ™ 4 . 
0, m+n _ odd 


Since (3) must hold for all A and arbitrary 6u;", the quantity in parentheses must 
i | & 


vanish and we arrive at the stress-equations of motion of order n 


° 
13, — Ts? +P = p Do AaMiu - 


m=O 
Strain. The three-dimensional components of strain, S;; , are expressed in terms of 
the u; by 
S;; = (u;,; + u;,;). (5) 


Upon substituting (1) in (5) we obtain, after a rearrangement of terms, 


co 


’ non) : 
Si; = Lai; ; (6) 
n=0 


where 
Si? = Bur) + us) + (nm + 1)(8,;0)""”? + 8;.05;"*"’)] 


in which 6,; is the Kronecker symbol. 
Stress-strain relations. In three dimensions, 


T = Ci 5x1 Sx ; c ji ™ Cia t= Ceri; > (7) 
S;; = 8i jx) kr , Sijkt = Sjikt = Skerii » (8) 
where c¢, ;,; and s;;,; are the elastic stiffnesses and compliances, respectively. 


The expressions for the two-dimensional 7’ in terms of the S{? 


inserting (6) in (7) and (7) in the first of (4), with the result 


Te 8.64: 2. Bakes (9) 


are obtained by 


Energy-densities. A strain-energy-density, U, and a kinetic energy-density, K, both 
reckoned per unit area of the plate, are defined by 


4 Y Y v(m y(n) 
U=h |] cinSe8er des = Feiss D> Dy Haw SF Si, 
A m=) n=0 
: m),(n) 
K =} | pu; U;,, a2 = tp >» >, Hat is 
h 


m=) n=U 


T = 0U/d8;; : 08;; /0S;; = 0, t et i. 
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3. Truncation of series. We retain only the zero and first order components of 
stress and strain and we write, tentatively, from (9), 


ry == Re. (10) 
Ti; = Fh'e, Si. (11) 


The components S‘° and S‘ involve only the displacements u{° , u{? and u{” and 
we neglect those of higher order. Then, following Cauchy, we neglect the velocity u3"? 
in the kinetic energy and equations of motion and provide for the free development of 
the strain S32? (= uj") by setting 732 = 0 in (10). The condition 73)? = 0 permits the 


22 
elimination of S;;’ from (10), with the result, again tentative, 
ery = 2hgiinrSee > (12) 
where 
Giixi = Ciser — Cij22 Coonr/C2222 « 


It may be observed that, because of the form of g, only five components of stress and 
five components of strain remain in (12). 

The first order terms, 7’! and S‘;’, are treated in a similar manner except that all 
three velocities u‘?) are neglected and free development of the three strains S}}’ is 
accommodated by setting 73!) = 0 in (11). Upon elimination of the S}} from (11), we 
arrive at 

4 Al ) 3 y(1) ‘ 
Ten Be YescaSet 3 a, b,c,d = 1,3, 
where 


tie, ol > | 
Yabca = A avea Sabed | 


in which |s,,.4| is the determinant 


| 
. 9 
Sia $1133. £8113 | 
' 
> i= . 2 9 
Sabcea | = $3311 $3333. £83313 


| 
| 


| 281311 281333 481313 | 
and A,,.4 is the cofactor of that element of s,,.4) which contains s,,. - 

As the final step in the process of truncation, the thickness-shear strains S};’ and 
SS. are replaced by k,S,)’ and k;S,;’ in the strain-energy-density, where k, and k; are 
correction-factors which may be used to adjust the thickness-shear frequencies to their 
exact values; thus compensating, in part, for the omission of terms of higher order in 
the series expansions. The correction-factors may be introduced conveniently by re- 
placing g,;,, with 


Qrinr = ke? j-okis1-29iinn (MO Sum), 
where m and n are the powers 

m = cos’ (ijr/2), n = cos (klr/2). 
Thus k”,,_, (or ki,,-2) is equal to k, , ky or 1 according as 7 + j (or k + 2) in gj jx: 18 
3, 5 or neither, respectively. 
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4, Recapitulation. The equations remaining, after truncation of the series and 
adjustment of the terms retained, are [7] 


Energy-de nsities*™ 


U = hgtiSi; Ser + Gh YercaSan Sea’, (13) 
K = ph(uprust + ahus us). (14) 
Stress-strain relations 
T° = 9U/aS\ = 2hg*nS, (15) 
TS) = OU/OS2 = HW yeeiSi (16) 
where 0S;;'/0S;; 0,2 ¥ j; ASS) /aSi? = 0, a # b. 


Stress-equations of motion 


7 ' r( . ( 
Tis ot F;” = Qhou;.e; 


; (17) 
Too ~ To + Fy” = Fh’ pus tr - (18) 
Strain-displacement relations 
Si = Bus) + ous) + 6203? + 8,05"), (19) 
Ss = uss + usd). (20) 
Note that, although Sj3’ = wu’) is contained in (19), neither of the two appears in any 


other of the Eqs. (13)—(20) because the zero order components of strain always occur 
as products with g*,,,; and the latter is zero when 7j or kl is 22. Hence, there are only 
eight components of strain and eight components of stress to be considered. The com- 
ponents of strain are related through the four compatibility equations 


Si 33 + $33.11 _ 2813.13 , S33 Tee Sis’s — S23)13 — Si2)33 ’ (21) 
S31 33 + S33 > 2833 13 Sits a Sst. = Sora o S32.11 ° 
Displacement-equations of motion 
Qhgtie(Ueers + Soxu;';) + Fj}? = 2hpuj’, , (22) 
2h Yaredls.aa — Zhghe(usrr + bu;") + Fy” = Fh’ pusre - (23) 


These equations are closely related to several predecessors. Thus, if the thickness-shear 
and flexure are eliminated by setting the transverse shear forces 73° and the couples 
and F;” equal to zero, the first and third of (22) reduce to the Cauchy-Voigt [1, 8] 


equations of low frequency extensional motion of thin plates: 


AN aredterag + Fs” = Bhpuy: - (24) 

Conversely, if the extensional deformation is suppressed by setting ui” = FS” = 0, 
(22) and (23) become equations of coupled thickness-shear and flexural vibrations [2] 

Dhgstese(Us sa + Uses) + Fo” = hous’. , (25) 

BN Varedtcaa — Lhgty(us, + us?) + Fy” = Bh’ puss: - (26) 


*Here and in the sequel: indices 7, j, k, 1 range over 1, 2, 3; indices a, b, c, d range over 1, 3. 
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Alternatively, in the case of isotropy, (22) and (23) become 
HUprea + (A + peace, + SRF? = purr , (27) 
k*w(uzrea + Uae) + 3h 2” = puree , (28). 
HUssca + (A! + p)ucien — Bh “hk wus’, + uy”) + FA OF? = pusir. , (29) 


where \’ = 2ud/(A + 2) and Xd and uw are Lamé’s constants. Equations (27)—(29) are 
also the isotropic forms of (24)—(26). Equations (27) are Poisson’s equations of low 
frequency extensional motion of thin plates [3 and 6, p. 497] or, what is the same, the 
equations of motion in generalized plane stress. Equations (28) and (29) are equations 
of flexural motion of isotropic plates with rotatory inertia and transverse shear defor- 
mation taken into account [9]. In the one-dimensional case they have the same form as 
the Timoshenko beam-equations [10] and in the case of equilibrium they have the same 
form as Reissner’s plate-equations [11]. Finally, by setting the transverse shear de- 
formation, S3°, and the rotatory inertia, 3h*pu ‘'), equal to zero, (22) and (23) may 
be reduced to equations equivalent to Cauchy’s [1]. 

5. Correction factors. The values of k, and k,; are found by equating the thickness- 
shear frequencies obtained from (23) with the corresponding ones obtained from the 
three-dimensional equations. 

In (23), let PF,” = us? = Oand u"’ = A exp iwt, where the A,” are constants. 
Then (23) become 

Ghi2aAg _ 4 ph" A,”. 


Upon equating to zero the determinant. of the coefficients of the A;"’, we have 


! | 2 3 
J2o20 — 9 Oa | = 0, g= t ph'w (30) 
which gives two frequencies, say w, and w, . These are to be equated to the two lowest 
roots of the frequency equation 
| C2321 o=.¢ 65, | = 0, c= 4ph’w/x’ (31) 


for the fundamental thickness-modes, as obtained from the three-dimensional equations 
[12]. The two equations yield a pair of equations which may be solved for k, and ky in 
terms of ratios of the ¢;;,; . 
As an example, consider the case of monoclinic symmetry with x, the digonal axis 
so that ¢2;22 and ¢2;.3 are zero. Then the roots of (30) are 
9: = Kigain = Keon , w, = 39,/ph’, 
2 2/ / ph? 
g3 = k3Qo323 - k3(C2303 — C2223; C2222) » 3s = 39s, ph 


and the two lowest roots of (31) are 


~) 


C, = Co: w, = rc,/4ph'’, 


/ 2 2 41/2 2 2 2 
Cz = 3{C2222 + Cos23 — [(Co222 — C2323) + 4C2223] “"}, ws; = WC3/4ph’. 


= x /12, k3 = ®'°C;/12go323 ° (32) 











56 R. D. MINDLIN (Vol. XIX, No. 1 
6. Straight-crested waves in a monoclinic plate. The adequacy of Eqs. (13)-(23), 
for the prediction of frequencies of vibration of bounded crystal plates, may be judged 
by a comparison of the five-branched frequency spectrum of an infinite plate, as obtained 
from (22) and (23), with the first five branches of the spectrum obtained from the three- 
dimensional equations. This is because the modes of vibration of a bounded plate are 
composed essentially of the anharmonic overtones of the modes of an infinite plate, as 
described in Secs. 16 and 17 of Ref [13]. Because of its importance in technology, the 
frequency spectrum of straight-crested waves traveling along the digonal axis in the 
AT cut of quartz is chosen as the basis of comparison. 

Quartz is a trigonal crystal (six elastic constants) and the AT cut is a plate which 
contains a digonal axis and whose normal makes an angle of 35°15’ with the trigonal 
axis [14]. When referred to rectangular axes in and normal to the plane of such a plate, 
the stress-strain relation has monoclinic symmetry (thirteen constants). With 2x, the 
digonal axis in the plane of the plate and x, the axis normal to the plate, the values of 
the thirteen constants, as computed from Bechmann’s [15] values of the six principal 
constants, are, in units of 10° newtons per square meter, 





Cin B64 C2303 = 38.6] C22 = —8.26 Coo33 = —1.42 
Cooon = 129.77 Cries, = 6651 Cus3 = 27.15 Coo03 = 95.70 
Caag3 = 102.84 Cyr = 29.01 Cys = —3.65 —Can03 = ~—«9.92 
C213 = 9.33 
while 
Cis = Co213 = Csais = Cosis = Cinn2 = C212 = Cazi2 = Coie = O. (33) 
Then, in (22) and (23), let F;” = F,’’ = Oand 
u;' = A;” exp i(fa, — wl), u,'’ = A,” exp i(éx, — wl). 


As a result, it is found that the first of (22) is coupled only with the second of (23) and 
the first (23) is coupled only with the second and third of (22). Thus, the general quintic, 
relating w and £’, reduces, in this case, to a cubic and a quadratic 
© 2 2 2 ° 
Jia — pw Oi s12k —19 foi 


oa 4 2 2 > — 
93128 Cisisgé — pw — 1g i312 ‘ie 0 


3th "Qeik 3h “Gisi2€ Yuanue + 3h "Seve = pw 


2 2 . aa 
Jigs — pw — 19h i236 -0 


3th “hi ask x + 3h "Gloss -~ pw 
The roots of these equations are plotted in ig. 1 to dimensionless coordinates 
Q = w/(x/2h)(ey212/p)"’’; @ = 2th/r 


The three branches marked flexure, face-shear and thickness-shear are the roots of the 
cubic, while the branches marked extension and thickness-twist are the roots of the 
quadratic. 

The five branches are to be compared with the first five branches of the analogous 
solution, by Ekstein [5], of the three-dimensional equations. The results of detailed 
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computations [16, 17] show correspondence to at least three significant figures over most 
of the range covered by Fig. 1. The asymptotic behaviors of the four branches which 
intersect at 2 =0,¢ = 0 


Q = $'(4/2)(y111;/Bei212) (flexure), 

Q = O(¥1313/¢1212) ” (face-shear) , 

Q = O(i111/C1212) (extension) , 

Q = —@'(4/2)¥1111/3ei212)'* ~~ (thickness-shear) 


are the same, for the approximate and exact equations, due to the method adopted in 
the derivation of the former. The remaining branch comes to zero frequency at 
ip = (Yii11€s en * = 0.7456 

in the approximate equations and 7¢ = 0.7467 in the exact equations. In the high fre- 
quency range the important frequencies, in the usual applications, are near 2 = 1 and 
the important branches are the thickness-shear and flexural branches. At 2 = 1, the 
thickness-shear branch is exact, due to the choice kf = 2°/12; the flexural branch in the 
approximate equations has @ = 1.2483 whereas the exact value is @ = 1.2417. 

The approximate equations should not be used for frequencies so high that branches 
higher than the fifth cannot. be neglected. In the AT cut of quartz the sixth branch has 
a real minimum at a value of Q slightly less than 2.0 [16]; however, due to the likelihood 
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of the contribution of an edge-mode from the complex segments of the sixth and seventh 
branches [17, 18], applications of the equations should be confined to 2 < 1.5. 


7. Uniqueness of solution. A theorem analogous to Neumann’s* leads to initial and 
boundary conditions sufficient to assure a unique solution of equations (15)—(20). 

Consider two systems of the five displacements, eight strains (omitting S32’) and 
eight stresses which satisfy the eight stress-strain equations (15)—(16), the five stress- 
equations of motion (17)—(18) and the eight strain-displacement equations (19)-(20). 
Let the differences between corresponding components of displacement, strain and 
stress constitute a ‘“‘difference-system” of those quantities. We form the equation 


| dt [ (7, + FO — 2hous,)u” 
to JA 
+ (TP, — TO 4 PO — BP pu?)u] dA =0 (34) 


in which all the components are those of the difference-system and t is an initial time. 
Now 


ghe [| at [ (uu + thu? uf?) dA = i) (K — K,) aA, (35) 
J to JA A 


where K is the kinetic energy-density of the difference-system and K,j is its value at ¢, 
Also 


Ts; uyer + (Tosa — Toe ust 
= (Tu) , + (Pet uf). — Te ul”, et T 2, a :, 
= (TP uf, + TP uf?),. — _ 1%(u + 62,u;"’),. be Us at 
= (T.; 7 0 + Tye Uest).a — ATS (Ui, + Us) + bo ,u; 
+ gin + 47S 3 + ui2)., , 
= (Ty 4 PVyD (8U/dS;) Si. — (BU/dSY) S32, , 
= oda + Tipu... —- Us; 


where U is the strain-energy-density of the liana Hence 
t 
ri 9) 0 7( ia) 
| at f (TP. + 72. — TH)us"] dd 
te JA 


- | at n (Tou + Ty) ds — | (U—U,)dA (36) 
to Cc 


JA 
in which the line integral is around the boundary, C, of the plate and the n, are the 
direction cosines of the outward normal in the plane of the plate. In the transformation 
from the surface integral to the line integral, it is assumed that the compatibility equa- 
tions (21) are satisfied. Finally, using (35) and (36) in (34), 


[ (U+K)dA = [ (Us + K.) dA + [ dt ll (Fu + Fu) ad 


JA 


e 


at 
o | ag nT ous + Tie us) ds. (37) 
/ to Cc 


*Reference [6], p. 176. 
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Then conditions sufficient for a unique solution are established by the usual argument 
based on the vanishing of the right hand side of (37) and the positive definiteness of 
U and K. It may be seen that, in addition to the initial values of uf”, uf”, uf? and u2? , 
there are five conditions to be specified at each point of the interior of a plate and 
five conditions at each point of the edge. In terms of orthogonal coordinates a, 8, 22 , 
the interior conditions are one member of each of the five products 

ca Se, Fa. Fa. FF 


In terms of orthogonal coordinates n, s, x. , the edge conditions are one member of each 


of the five products 

ry (it ) ry’ ks ry (0 (0 y(1). (1) =f ) 

“a. pa ae ee ae i. . ea (38) 
The first two of these products give the boundary conditions of generalized plane stress 
and the last three give the boundary conditions of Reissner’s equations of flexure [11]. 


8. Orthogonal functions. A theorem analogous to Clebsch’s* may also be established. 
Consider two solutions 


(us? us’) = (u;”", ws’) exp ta, f 


i 


I] 


(0) ) 0) (1) . 
(us, uy’) = (us?", uy") exp iw, t 


of the homogeneous (F;" = F;"”) = 0) stress-equations of motion, so that the equations 
2phau;?” = —TiH i, 
2 ohrwius?” = —TY + Tr” 
—2phou;? = Ti); ; 


—3ph*wiu,?’ = TS’ — Tr 


(O)e (l)e (O)r (1)r 


are satisfied. Multiplying these equations by u;"’*, u,'’", u;"’", u,’”, respectively, adding 
and integrating over the area of the plate, we sheiin. on the ‘eft aad side, 


2hp(w _ w*) (u\° “uc? . 4. 1p” “— a) dA 
JA 
and, on the right hand side, 
rr * O)r my (O)r O)s ryvil)e l)r pie,, a (0) r l)e WSG7 0, SE) 
| (T; ;u — 7 ,u +. 7 <. alls = To al a 7; an Us aad T's) Uy, ) dA. 


The latter, by a process similar to that employed in the preceding section, may be 


transformed te 
2h | [gh SP Si" — SHPSIP*) + AW Yerca(Sar*Sia” — SSP") dA 
4. g n (Ti? us Boon. per 0) ” 0° . “4. oo, 8 cae po l)r <) ds. 
. 
The integrand in the surface integral vanishes identically. The integrand in the line 


*Reference [6], p. 180. 








60 R. D. MINDLIN [Vol. XIX, No. 1 


integral vanishes for any one of the thirty-two sets of edge-conditions obtained by 
equating to zero one member of each of the five products in (38), or for any of the associ- 
ated conditions of elastic support. Then, if w, ¥ w, 
| (uj; "uj; * + thu,’ ’u,'’*) dA = 0. 
JA 
Some observations may be made regarding normal modes of rectangular plates. In 
the general triclinic case there are no solutions of the displacement-equations of motion, 
(22)—(23), with each of the five components of displacement a single product-function. 
However, such solutions are possible for monoclinic and higher symmetries. 
Thus, in the monoclinic case, where (33) hold, (22) and (23) are satisfied by the 


four sets of displacements (and no others) 


u, = A, cos (&, + pr/2) sin (fa; + qw/2) exp iwt, 
us = A,” sin (fx, + pr/2) cos (fx; + gr/2) exp iwt, 
us. = Ags’ sin (fx, + pr/2) cos (fx; + qr/2) exp iwt, (39) 
u;'’ = A," cos (fa, + pr/2) cos (fx; + gr/2) exp iwt, 
u;'’ = A," sin (€, + pr/2) sin (tx; + gr/2) exp iwt, 


p = 0,1; q= 0,1 


subject to a quintic equation relating w, & and ¢°. Fora plate z, = + h,,x7, = +h;, 
normal modes, in which each displacement is a single product function, are possible if 
and only if, onv, = +h,, 
0 (0 l mio ryyviii 
UW =u =u =Ty, =Ty, =0 
(40) 


or T,. =T; =T;; = “= ney = (0 


and, onxz; = +h;, 


l =u =u, =T7Ti; =T33 = 0 (41) 


whence 


£ = mrx/2Qh, , ¢ = nw/2h; ; mn=0,1,2--- (42) 


and the frequencies are given by the five roots of the quintic for each m and n. 

The first set of (40) and the first set of (41) are the conditions analogous to simply 
supported edges in the elementary theory of flexure; but note that, in the present case, 
(41) are not obtained from (40) by interchange of indices 1 and 3. As in the elementary 
theory, solutions in closed form are not possible, in general, if all four edges are free; 
but closed solutions may be obtained with one pair of parallel edges free and the other 
pair under conditions (40) or (41). Then each of the displacements in (39) is given by a 
sum of five functions and one of the sets of roots (42) is replaced by the roots of a5 X 5 
transcendental, determinantal equation. There is one exceptional case of a solution in 
closed form, with all four edges free [19], analogous to the Lamé modes of an isotropic 


plate. 
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(Continued from p. 50) 


The author does not discuss these matters at all. This is unfortunate, since there is little value to 
steady-state analysis unless one shows that the dynamic process asymptotically approaches the steady- 
state process as the length of the processes increases. Furthermore, it is essential to indicate the rate 
of approach. 

The author sets himself the task of determining steady-state policies under the assumption of their 
existence. Granted the existence of a ‘‘steady-state,’’ the function f;(n) has the asymptotic form ne + 
b; + 0(1) asm — ~ where c is independent of 7. The recurrence relations then yield a system of equations 
for c and the };. 

This system can be studied by means of linear programming as a number of authors have realized; 
cf. A. S. Manne, ‘‘Linear programming and sequential decisions,’’ Management Science, vol. 6, 1960, 
pp. 259-268. 

Howard uses a different technique based upon the method of successive approximations, 
an approximation in policy space. It is a very effective technique as the author shows by means of a 
number of interesting examples drawn from questions of the routing of taxicabs, the auto replacement 


in this case 


problem and the managing of a baseball team. 
The book is well written and attractively printed. It is heartily recommended for anyone interested 
in the fields of operations research, mathematical economics, or in the mathematical theory of Markov 


processes. 
RicHarD BELLMAN 


Handbook of supersonic aerodynamics—shock tubes. By I. I. Glass and J. Gordon Hall. 
Section 18. Supt. of Documents, Washington 25, D. C. xxxviii + 604 pp. $3.75. 


An extensive handbook of shock tube technique. The first part, on the theory and performance 
of simple shock tubes, is by I. I. Glass. The second part, on the production of strong shocks and the 


application, design and instrumentation of shock tubes, is by J. G. Hall. 
R. E. Meyer 


The theory of thin elastic shells. Edited by W. Koiter. Proceedings of I. U. T. A. M. 

Symposium, Delft, 1959. North-Holland Publishing Company, Amsterdam, 1960. 

ix + 496 pp. $9.00. 

The volume contains the papers presented at the Symposium on the Theory of Thin Elastic Shells 
sponsored by the International Union of Theoretical and Applied Mechanics (I1.U.T.A.M.) and held 
at the Technological University of Delft in August 1959. All papers at the Symposium were delivered 
by invitation only and attention was confined to two particular aspects of shell theory, viz. nonlinear 
theory and problems lacking axial symmetry. An Introduction to the Symposium by Prof. C. B. 
Biezeno is followed by the contributions of the following authors: W. T. Koiter, W. Zerna, D. G. Ashwell, 
B. Budiansky, H. Ebner, E. I. Grigolyuk, J. H. Haywood, G. Czerwenka, W. Schnell, M. Kuranishi 
and J. Niisawa, A. L. Bouna, Y. Tsuboi and K. Akino, A. van der Neut, I. N. Vekua, M. W. Johnson, 
P. M. Naghdi, W. A. Nash and J. R. Modeer, B. R. Seth, P. Seide, N. J. Hoff and J. Singer, J. W. Choen, 


E. Reissner, H. Neuber 
R. T. Sareip 


SOLUTION OF A DIFFERENCE EQUATION PERTINENT TO LINEAR, 
PARAMETRIC ELECTRIC NETWORKS* 


By B. J. LEON (Hughes Research Laboratories Malibu, California) 


This paper presents an algorithm for finding the fundamental system of solutions of 
a class of linear difference equations relative to a theorem of Perron [1]. This particular 
fundamental system is important in the analysis of physical systems with periodic 
parameters such as linear parametric networks [2] and parametric amplifiers. These 
difference equations arise as the Laplace transforms of the linear differential equations 
with periodic coefficients that characterize parametric networks. These transforms must 
be analytic in a half plane. In Ref. [2] it is shown that the method of variation of param- 
eters yields such an analytic transform if the solutions to the homogeneous equation 
are Perron’s fundamental system. These solutions to the homogeneous equations need 


not be analytic or even continuous. 
The homogeneous linear difference equation under consideration is 


0 = Uwtn) +4,-,.@Ue +n — 1) +--+ + a,@WU@) (1) 
subject to the conditions 


a) lim a,jw+m) = A;, a constant, 


m—-@ 


integer 


b) A, #0, and a(w+m) #0 for m=0,1,2,-:-, L (2) 
c) u" + Ayu" + ess + Ay = (Cu — u)(CU — UW) +++ CU — 4,) 





d) |u;|lu;| for tj. J 


In the difference equation literature—Milne-Thomson [3] is the most complete 
reference for methods of solution—the emphasis is on finding analytic solutions to the 
homogeneous equation. The complexity of these methods increases markedly with both 
the order of the difference equation and the complexity of the coefficients of the equation. 
With these methods, only the simplest of parametric amplifier problems can be put on 
an IBM 709 computer. If a few second-order effects are included in the amplifier circuit 
model, the machine overflows. 

The algorithm presented in this paper is intended for execution on a digital computer. 
When the problem is set up on a digital computer, the coefficients in Eq. (1) are sub- 
routines, which are brought into the main program as numbers. The number of operations 
in the main program depends on the order of the difference equation and on the asymp- 
totic behavior of the coefficients. The complexity of the coefficients for small w does 


not affect the complexity of the main program. 


*Received August 12, 1960. 
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The basis for our algorithm is the property of the solutions to the difference equation 
given by Perron’s theorem. That theorem states that for Eq. (1) subject to the conditions 
(2) there exists a fundamental system of solutions, U, , U,, --- , U, , such that 
Uw + m + 1) 


lim — 
ane U (@ + m) 


integer 


= u 


. (3) 


We recall that the definition of a fundamental system of solutions for difference equa- 
tions is similar to the definition for differential equations except that the arbitrary 
constants in the differential equation case are replaced by arbitrary periodic functions 
of period one in Re [w]. 

The algorithm has five basic computational steps that are repeated over and over 
until all the desired solutions are found. Let us first describe the steps in words and 
then present them formally with justification. The steps are 

(I) There is a straightforward procedure for constructing solutions to Eq. (1) [4]. 
Although the solutions so constructed do not have the asymptotic property (3), we may 
formally write each solution as a linear combination of the members of the desired 
fundamental system. 

(II) From Eq. (3) we see that the asymptotic properties of the desired U’, are all 
different. Thus if we examine solutions generated in step (1) for large values of the 
argument we can ascertain their relative dependence on the desired U; which has the 
fastest rate of growth. Call this largest solution U, . 

(III) Starting from n-linearly independent solutions constructed in step (I) with 
known U, dependence from step (II) we can proceed by systematic elimination to find 
(n — 1) solutions which depend only on U,, U;,--- ,U,. 

Since all the | uw, | are distinct, steps (II) and (III) can be repeated eliminating the 
U, dependence, then U,; , etc., until there remains only one solution which has the 
asymptotic property of U,, . (We choose the subscripts so that | u, | > |u| > --- > 


ts |) 


(IV) When a solution to an nth-order linear difference equation is known, that 
solution can be used in conjunction with the equation to derive a new linear difference 


equation of (n — 1)th order [5]. We shall show that if the solution U,, , found in (III), is 
used to derive an equation of order (n — 1) then this equation has a fundamental system 
of solutions UU; _, with 
‘ (wo + m + 1) u 
1\W i ‘ 
lim — =—, i = 1, 2, ,(n — 1) 
igae Us + m) Uy 


The U;,, , like the U; , have distinct asymptotic properties. Thus we can proceed by 
steps (I), (II), and (III) to find a solution U,_,,, . This solution can then be used to find 
an equation of order (n — 2) and so on. 

(V) Finally we get the desired U,; from the solution U;,, by 


Uw) = U,() S: [U,..(X)] AX, 
where S“[ | AX is Norlund’s summing operator [6]. The various steps of the algorithm 


are repeated until all the desired U; are found. 
Formally the algorithm may be written as follows: 
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I) Let 
a for 0 < Re [a] <n 
| —fa,i(w — n)U'(w — 1) + a, — n)U'@ — 2) + --- 
| + alw — n)U'w — n)| recursively first 
for n < Re [fw] < (n + 1) 
then 
U'(w | (n + 1) < Re [w] < (n + 2) 


ete. for all 


Re [w] > n 


[U'@t+tn) +a,i@U'w+tn—-—1)4+-::: 
+ a,(w)U ( + 1] recursively first 
for 0 > Re fw] > -1 
then —1 > Re [fw] > —2 
etc. forall Re [fw] < 0. 
Justification. 


A) The U” are well defined solutions by condition (2b) and Ref. [5]. 
B) The (” form a fundamental system of solutions because 


a) for0 < Re [w| < 1, Casorati’s determinant, 
I w wa”! 
l w l (w i. 
Dw) = + ? 
lL wtn-1l «+ wWwtn-1)"' | 


is a Vandermond determinant [7], which is easily seen to be non zero. 
b) By Heymann’s theorem [8] 


Dw + 1) = (—1)"a) D@). 


c) By Condition (2b) above D(w) # 0 for all w. 
d) By Casorati’s theorem" 


D(w) # 0 = the U's form a fundamental system. 
C) The U’’s can be written formally as 


vu’ _ pU, + p.U, +: + DiUn ’ 
U? = p'[pU, + pie + +++ + piUn), 
U* = p'[pU, + p.U2 + +++ + p.U,], 


U" = p"[pU, + pU. + --- + piU,l, 
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where the p’’s and p;’s are periodic in Re [w] with period one. The u,’s are ordered so that 


it ite tf aoe > 18 


' = 3? 


(II) If U'\(w + m) ¥ O let 





‘w, m) = Uwt m) 
thst di U*(a + m) 
Justification. 

If p(w) + 0 then 


lim p'(w, m) = p(w). 


To show this we note 
a) By (3) for every « > 03 .M;(w) and a real, positive K;(w) 3 form < M,(e) 





, \/! \ m— Mi (w) | 7 , \/| m—M¢ 
K,(o)(| u; | — 0"? < | Ui@+ m) | < Ki@)(|u; | + 9. 
b) For 
[uy | — | ue | 
e< 5) 
and m > max J; M, we have 
| Uw + n) Kw) (ju; | + eed | u;|+e sutton 
Tri ee eae < - ss | M VW | rg yee: : 
U,(w + m) | K,@) (| u, | — eer u,|—e 
So for] ¥ 1 and K,(w) ¥ 0 we have 
lin Uo + m) 
m ———_—_——— = 
mao Ui(w + mM) 
c) The rest is straightforward. 
(III) Let 
1 U' 1 . 
=—-—l for += 2,3 n 
p* b 
Justification. 


The U*"’ can be written formally as 


ee p; U2 + Ds 'Us i D, gm ’ 
1] 2,177 


p ‘[p2 /2 + p,'Us; aa eee + p."U,); 


~~ 
i 
eS 


U"* = p™"[p2"U2 + p3'Us + «+: + pr'U,), 
where the p''’’s and p}'’’s are rational combinations of the p'’s and p}’s. 


This form is the same as I(c), so steps (II) and (III) can be repeated until there 
remains only 


‘Ogee coo D, ues im : 


This solution is a member of Perron’s fundamental system relative to the root u, 
Henceforth we suppress the superscripts and refer to this solution as U, . 
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Consider the equation 


U s(w +n -— 1) + a, 2,1(w) U sw +n 2) +--+ + i(w) Uw) = 0 (4) 
with 
— A, ,1\w) 


= a(w)U,(w) + --- +.a,(w)U,@ + 8). 


Justification. 


\) Equation (4) has a fundamental system of solutions [5] 
i ae t#=1,2,-- ,m—1})), 
U, 
whe re 
AMX) = f(X + 1) — f(X). 
» tee tete Ss 


U, (we + m) U,, 


To show this we write 


I w+ m+ 1) E (w + m + 2) bay U,(w ttt 
U’. w + m U.(a + m + 2) Ufo + m + 1) 


Uo +m+i1)_ U,@ + m) 
{| Ze +m+ 2) /| eet m + 0 
U (ew + m) U,@ + m) 
_ | Hale +m + 2) / | wet m+ 
U,(@ + m) U,(@@ + m) J 


{([ ue +m +0} /| Ue +m +D)) 7 i} 
( U (a + m) U,(@ + m) 


| Uae t m + 1) U@t “| 


——.. 


lin U. (oa + m+ 1) (u?/u2) — (u;/u,) u, 
im - - "=> . 


rages U; ,(@ + m) Pr (u;/t,) — 1 Me 
C) The construction procedure (I) gives a fundamental system of solutions to Eq. 
1). By Heymann’s theorem the only zeros of D(w) are congruent to the zeros of ao, , 
but these are considered singular points of the equation. Thus the hypotheses of Casorati’s 
theorem are satisfied. 

D) Repeated application of (II) and (III) leads to a solution with the asymptotic 
property of U,,-:)., . This solution, along with Eq. (4), can be used to construct an 
equation of order (n — 2) with the same properties as Eq. (4). 
V) Let 


U; = U, S? [U;.,(X)] AX, 


where {2 [{ | AX is Norlund’s summing operator [6]. 
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Justification. 

These U; are the solutions U; to Eq. (1) that we desire because {* 
A (step (IV)A) to within a periodic function. Thus the problem is solved. 

Discussion. In the various steps of the algorithm divisions are required. For these 
operations we must be sure that the divisor is non zero. First we note that if the co- 
efficient of Eq. (1) are analytic along a line parallel to the real w axis except for isolated 
singularities, then the solutions U' generated in step (I) will have the same property 
along that line except at the integer values of Re [w]. Therefore the zeros of the various 
divisors will be discrete, and the algorithm can fail because of attempted division by 


is the inverse of 


zero only at a discrete set of points. 

The other place where the algorithm as stated may fail is in step (II) where we 
require pi(w) * 0. It is possible that our construction of the U' gives a U* which is 
independent of U, . Since the U' form a fundamental system, at least one of the UL” 
must depend on U, . Therefore if we find that the p‘(w, n) do not tend to a limit, we 
must renumber the U’' so that the new U’ is one of the U‘ that depends on U, . Through- 
out the elimination process there is the possibility of this analogous situation occurring. 
However, it can always be circumvented by renumbering the U’. 

In the parametric network problem there are symmetries in the original homogeneous 
equation. These symmetries can be used to reduce the number of operations required 
in the algorithm. In a subsequent paper these symmetries and the corresponding re- 


ductions will be discussed. 
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A NOTE ON THE MINIMAL REPRESENTATION 
OF TRIGGERING MATRICES* 


BY 
ARTHUR GILL 


(Department of Electrical Engineering, University of California, Berkeley, Calif.) 


Abstract. A triggering matrix is an array of numbers in which each row represents 
an actuating input of a data-processing system. In this note, upper and lower bounds 
are given to the minimal form of such a matrix. 

In this note the term ‘“‘triggering set’’ will refer to any group of differing numbers, 
each of which serves as an input to a computing or control system. The members of the 
set, called ‘‘trigger-words”, have no numerical significance—their only function being 
the initiation of various processes in the system. It will be assumed that each trigger- 
word initiates a different process, and hence that the only requirement for satisfactory 
system operation is the ability to differentiate among the words of the set. Typically, a 
triggering set may be the digital representation of analog quantities to be fed into a 
recognition system. 

In the following discussion, m will denote the number of trigger-words in the set, 
n’ the number of digits in each word, and r the ‘‘set radix’’—the number of values that 
each digit can assume. m, n’ and r will be taken as finite. For convenience, a triggering 
set will be written as a matrix, called the “triggering matrix’”’, in which the element 
common to the 7th row and jth column will represent the jth digit in the 7th word. 
Matrix (1) shown below is an example of a triggering matrix with m = 6, n’ = 9 and 
} 3. 


123456789 
212011201; 


2)210110220) 
3,;102220001 


~ 


| (1) 
4;100211121] 


5'02110211 2! 
61021022110) 


A problem which often arises in data-processing systems with inputs characterizable 
by triggering matrices, is that of determining the smallest number of digits necessary 


*Received November 16, 1959; revised manuscript received January 8, 1960. 
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for satisfactory system operation. In terms of the triggering matrix, it is desired to find 
the smallest number of columns such that the rows are still distinct. Under noisy con- 
ditions it is further desired to maintain a specified minimal distance between the rows 
in order to increase reliability (“distance”’ having the usual Hamming connotation, 
extended to codes of radix r). Deleting the largest number of columns from a triggering 
matrix such that the distance between any two rows is not less than d, results in a matrix 
which will be called the “minimal triggering matrix consistent with d’’. In a previous 
paper [1] the bounds on this deletion were produced for the special case r = 2 (i.e., for 
binary triggering matrices). The following is the extension of the previous results to the 
general case. 

Theorem. Let the dimension of a triggering matrix of radix r and minimal distance 
d, be m X n’. Let the dimension of the corresponding minimal triggering matrix con- 
sistent with the minimal distance d < dy be m X n. Then 


E =1) + jlog, m} <n < d(m — JI), 


where [x] denotes the largest integer which is not larger than x, and {x} the smallest 
integer which is not smaller than ~. 

The lower bound equals the minimum number of columns which can accommodate 
radix-r rows with minimal distance d. This number can be found by generalizing relation- 
ships previously derived in [2] and [3] for the binary minimal distance codes. For that 
purpose, let a be a sequence of n digits of radix r 


a = (a, , Qe, °*** Qa) a, =0, or 1,°---, o r—l. 
The sequence a © £ will be defined by: 


! 


aOB = (|a, — B|,|a2 — Bl, +--+, lan. — Ba | 


The distance between a and 8, denoted | a © 8 |, is then: 
jJeOB| = Lila: — BI. 


Now, let J/ be an m-word n-digit code of radix r and minimal distance d. Let 


a, B (a # 8) be any n-digit words of radix r which end with n — [d — 1/r — 1] zeros. 
Then: 
' d—1 
aati <¥ «0h —41 ed ~ 4. 
r—1 
If y, 6 (y # 46) are any elements of M, then y © a and 6 © 8 cannot be identical. Conse- 
quently: 
mept rm < ® 
or: 


n> E — | + log, m. 
r—l1 


Since n must be an integer, the lower bound of the theorem follows. 
To verify the upper bound, consider an m X n’ {log.r} binary matrix produced by 
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replacing each element in the original triggering matrix with its binary equivalent. 
This is done in matrix (2) for the ternary triggering matrix (1). 


ye 


' 2 3 Se @& @ @ 8 
(TOOT (565661 67 166067] 
100100010100101000 (2) 
010010101000000001 
4/01000010010101 1001 
Ss |/001001010010010110 
« [(001001001010010100, 


oe = 
° 
° 








In the new matrix there are n’ groups of {log,r} columns, each column group repre- 
senting a single column in the original matrix. Now, from previous results (see [1]), it 
is known that the minimal binary matrix consistent with d contains at most d(m — 1) 
columns, and hence at most d(m — 1) column groups. Converting these selected column 
groups back to their original radix-r representation results, then, in a matrix of radix r 
which contains not more than d(m — 1) columns and whose minimal distance is at least 
d. The upper bound is thus verified. 

Both the lower and upper bounds in the theorem are achievable with equality. The 
lower bound is attained, for instance, when m is a power of r, d = 1, and the triggering 
set contains all possible radix-r words of length log,m. A triggering matrix in which the 
upper bound is attained is shown in (8). 








o— m-| —oe— m-i— —nm-i— 

1oO- ooio.- OO+rrteer 1o ° 

O1r-+-oo00! OOreeeee Oreo 
Mile ccceereseeceesssesecese 

O00++1000++tOrereeee8 00++10 (3) 
O0+*O100++Olereee 00--O! 

O0+--0000 OGercsecs 

é(m-1) 


It can be concluded that both bounds on the minimized number of digits are inde- 
pendent of the original number of digits, and that the upper bound is also independent 
of the number of values which each digit can assume. 
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ON OPTIMUM RECTANGULAR COOLING FINS* 


BY 
CHEN-YA LIU 
Carnegie Institute of Technology Pittsburgh, Pa. 


1. Introduction. In the past, most discussions concerning the cooling fins of a heat 
exchanger were confined to cases where the heat from the fin surface was eventually 
dissipated by the surrounding convective fluid. For this case, the governing equation 
is linear and the solution can be obtained without difficulty [1]. However, as space 
technology advances, a heat exchanger may have to be designed for an environment 
where the only heat transfer mechanism is by radiation. Furthermore, in any space 
vehicle design, the over-all weight of the vehicle is of utmost importance. It is, therefore, 
desirable to know the fin geometry of least weight. The essential difficulty in dealing 
with cooling fins when the convective transfer mechanism becomes insignificant arises 
from the fact that the governing equation is no longer linear. Up to the present, the 
solutions for this type of problem were obtained mainly by numerical method [2, 3]. 
It is difficult to optimize the result thus obtained. 

The present paper presents a parametric solution for a rectangular cooling fin in 
terms of known functions, from which the optimum geometry of the fin with least 
weight is uniquely determined. 

2. Statement ofthe problem. The governing differential equation of the temperature 
T(x) in a rectangular thin fin as shown in Fig. 1 is 


aT ain 
a (1) 




















Tl, 1 b 
L 
Fia. 1 
where 
s = heat transfer coefficient, 
k = conductivity of the fin material, 


b = fin width, 
= a constant, equal to | — 4 in actual application and assumed to be greater than 


unity in the present problem. 


*Received March 8, 1960. 


1961) NOTES 73 


The boundary conditions for a long fin are 


Pf, & o06 (2) 
aT 
7 0 at s=L. (3) 
The cooling rate of the fin in terms of the rate of heat conducted out of the fin base is 
aT 
q —_ —kb dz Lae (4) 


The problem may now be stated as: For a fixed fin weight (or area A = bL), determine 
the value of b such that g is a maximum subject to the governing equation (1) and the 
boundary conditions (2) and (8). 

3. The temperature distribution. Though the solution of the system of Eqs. (1), 
(2) and (3) can not be obtained in closed form, it is possible to derive a parametric 
relation between 7 and x as shown below. 

Multiplying Eq. (1) by d7'/dx and integrating gives 





ary’ — = 
(= “cin **"* (5) 
The constant C is determined by Eq. (3) and by assuming 
T=T, at z=L. (6) 
Then, 
Re eee, «a e 4s na +1 rd 
°-Cim * @ 


Note that the boundary value problem is completely defined by Eqs. (1), (2) and (3). 
The additional condition, Eq. (6), is used merely as a parameter which shall be uniquely 


determined. From Eqs. (5) and (7) 


aT 7 , 4s 1 ror (mp 
dx lz + “5 | : : T : (8) 


Introducing the new variables 





t = (T,/T)*", (9) 
to = (T1/T»)*"". (10) 
Equation (8) becomes 
/ 1/2 
?(1 — 2) va a - | ete ee" (11) 


where 8B = (a — 1)/2(a + 1). 
Integrating Eq. (11) gives 


a B-1/ —1/2 4( + 1)s gee / 7 d 
[ Pa- 1 a = [sets] ris? | de (12) 
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or 


B(8, 2) — B.(6, 2) = E + Ne] Te-MAL, — 9), (13) 


where B is the complete Beta function and B, , the incomplete Beta function which is 
defined for 0 < t < landg>0O [4]. The two conditions for the existence of the Beta 
function will be examined with regard to Eq. (13): 

(i) when ¢ = 0, Eq. (13) fails to give a meaningful solution. However, for this case, 
either 7’ — « which is physically impossible or 7’, = 0 which will be discussed at the 
last section as a special case. For the present discussion, we assume 0 < ¢ < 1 in Eq. (13). 

(ii) when 8 > 0, a > 1. If a = 1, the problem is for the case with the ordinary 
convective surface condition, for which the solution is available [1]. 

Under the restrictions of 0 < ¢ < 1 and a > 1, Eq. (13) expresses the functional 
relation between 7’ and x with 7’, as the parameter. Furthermore, the values of T and x 
have a one-to-one correspondence, since both sides of Eq. (13) are single-valued functions. 
To determine the parameter 7’, , Eq. (2) is substituted into Eq. (13). 


’ s i 1/2 yes : 
B(6, 3) — B,.(8, 3) = | 7 ee | Tia-v?, (14) 


After 7’, is calculated from Eq. (14) for a fixed physical condition, the rate of heat 


transfer, Eq. (4), becomes 


a Fe eine 
gai —— ft °° ~T. ) (15) 
f FE + 1) . 
4. The optimum geometry. Equation (15) gives the rate of heat transfer in terms 
of the parameter 7’, which has to be determined from Eq. (14). For a fixed A = bL, 


T,, will vary with b, so that q is, in general, a function of b and 7, . If Eqs. (14) and 
(15) are rewritten as 


q(b, Tr) = [ | ra. << S) (16) 
\a + 1 


-.. pear 

p(b, Tr) = B(B, 4) — B,.(8, 3) — [ see) rere? = 0, (17) 

where q and p are now to be considered as functions of b and 7’, , the equivalent problem 

of optimizing q will be to find the extreme values of the function q(b, T',) subject to the 

subsidiary condition p(b, 7;) = 0. The stationary point thus found gives the desired 

optimum values of b and 7, . The solution to this problem can be obtained by means 
of Lagrange’s multiplier m defined by [5] 


oq Op ai 
2 — 8S 
0b sas ob 0, om) 
oq Op 
— —A- == (), (19 
OT , sae OT’, } 


Equations (17), (18) and (19) serve to determine the stationary values of b, 7, and the 
constant m. Eliminating m from Eqs. (18) and (19) gives 





~I 


or 
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bY? -[¢ 7 vee | [15 "' = 22a + Ie = TE", (20) 
(a + DkTo Gy" + TF5"" 


Comparing Eqs. (16) and (20), it is seen that tf) > (a — 1)/2(2a@ + 1). Though Eqs. 
17) and (20) are sufficient to determine the optimum value of } by eliminating 7, , 


it is more convenient to write Eq. (20) as 


lide — 1) —G 4G {24a + 1) + GY), (21) 


S(2a + 1) 


\ here 
G = (a + I)kb*/(2a + 1)sA°TS™' (22) 


and the positive sign in front of G’ , is chosen because i, >(a- 1) /2(2a + 1). Equation 


17 becomes 


B(8, 3) — B,,(8, 3) = 2@ + IQa + 1)7'°a"6 , (23) 


where /, is given by Eq. (21). It is seen that the transcendental Equation (23) contains 
only the variables a and G. For a fixed a > 1 and under the conditions G > 0 and 4, > 
(a — 1)/2(2a@ + 1), Eq. (23) is satisfied by one, and only one, value of G which, in 
turn, determines uniquely the optimum geometry of the fin from Eq. (22). 

It remains to be shown that the extreme value of the function q(b, 7.) does exist at 
the stationary point thus found and that the extreme value is a true maximum. For 
the existence of the extreme value, it is necessary that the two partial derivatives dp/db 
and @p/dT, shall not both vanish at the stationary point [5]. It can easily be shown 
that both derivatives will vanish only when 7, = 0 which is excluded in the present 
problem. The sufficient condition for the extreme value of q to be a maximum is that 
dq < 0 at the stationary point. This is indeed true if tf) > (a — 1)/2(2a + 1). 

5. Special case when 7', = 0. In the above discussion, the case when T, = 0 
is excluded. If 7, = 0,C = 0 from Eq. (7). A direct integration of Eq. (5) with boundary 


condition (2) gives 


my—(a—-1)/2 fete (a — 1)’s ac 
, i = [‘ + us| . (24) 


— a ae 2 
q = (- +) 7 ; (25) 


It is seen that no maximum value of q exists for any finite value of b. The non-existence 
of the extreme value is apparent since physically 7’, = 0 requires the fin length L — «. 
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GENERALIZED THERMAL RESISTANCE 
BY 
J. A. LEWIS* 
Bell Telephone Laboratories, Inc. Murray Hill, New Jersey 


1. Introduction. The thermal resistance of a conducting body V is usually defined 
(see, e.g., [3]**) in terms of a harmonic function u(2, , x , x3), satisfying the boundary 
conditions 


u,, on 8, 
u=4 


lo, on So, 
un, =0, on S,, 
where S = S, + S, + S, is the boundary of V, n; is the outward normal vector, and 
u, and uy are constants. We shall refer to S, as the inlet, S, as the outlet, and S, as 


the adiabatic surface of the body V. 
The thermal resistance R of this configuration, i.e., the body V with given inlet and 


outlet, is then defined by the equation 
R(u) = (u, — up)/Q(u), (1.1) 


where 


Qu) = | un; aS, 


the total transmitted heat power. 
By applying Green’s identity to Eq. (1.1), using the harmonicity of u in V, we obtain 


the alternate forms 


Rtu) = (uy — u)°/D(u) = D(q;)/Q°(4q:), (1.2) 
where gq; = u,; , Q(q;) Q(u), and the Dirichlet integrals D are given by 
Du) = D(q) = | agi dV = | uu, dV. 
JV V 
If %,; — Up is interpreted as a potential difference and Q as total current, Eqs. (1.1) 


and (1.2) also give the electrical resistance of the configuration for unit electrical con- 
ductivity. The two forms given by Eq. (1.2) are simply the familiar “2°/R’’, “7°>R” 
relations between the power dissipation D, the potential drop FE = u, — uw, , and the 
total current, J = Q(u). 

Since u; — wp is obviously the maximum temperature difference in the body V, R(x) 
gives the maximum temperature difference for unit transmitted heat power and thus 
serves as a power rating for the body, at least for these boundary conditions. Unfortu- 


*Received April 2, 1960. 
**Numbers in square brackets refer to the bibliography. 
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nately, in contrast to the problem of electrical conduction, in practice one encounters 
a wide variety of boundary conditions in heat flow problems. We shall show, in fact, 
that the temperature rise in general is greater for given input heat power than that 
given by the above Dirichlet conditions, so that R(u) is unsafe to use as a power rating, 
unless, of course, it is certain that the above boundary conditions are satisfied. 

Very often in practical problems the exact boundary conditions are unknown. In 
the following, we consider two other boundary value problems, which we regard as 
typical, namely the Neumann problem and the Robin-Neumann problem, both with 
piecewise constant boundary values. We extend the definition of thermal resistance to 
these cases in a natural and expeditious way and obtain upper and lower bounds, which 
include the resistance in all three problems between them and thus serve as universal 
bounds, at least for the class of problems considered. 

2. The Neumann problem. The temperature v, again harmonic in V, now satisfies 


the boundary conditions 


a ’ on S, ’ 
yn; = e on S,, 
— Q/So ’ on So ; 


where Q is the total input heat power. In this case the thermal resistance is not usually 
defined. We may extend the definition, however, by using the analogy with electrical 
conduction and requiring that the resistance R(v) again satisfy the usual “J’R” relation 
between the power dissipation D and the total current Q, namely 


Dive) = QR). 
If we apply Green’s identity, we obtain 
Rv) _ ((v), = (v)o)/Q, (2.1) 


where the average inlet and outlet temperatures (v), and (v)) are given by 


aa 8," | v dS, (v)o _ so* | v dS. 
JS Se 


As before, we have the alternate forms 


Rv) = (v), — @&)o)’/DO) = D(q:)/Q, (2.2) 


where now q a 

If we compare Eqs. (1.1) and (2.1), we see that we can include both in a single 
definition, if we define the resistance in general as the ratio of the average inlet-outlet 
temperature difference to the total input heat power. This is the definition which we 
shall adopt. It is particularly apt because it admits immediately of a simple physical 
interpretation and because, at least in the two cases considered so far, it may be ex- 
pressed in terms of the Dirichlet integral which is easily bounded. 

3. The Robin-Neumann problem. We now consider a heat flow problem which we 
regard as canonical and from which the two preceding problems may be derived as 
special cases. We again seek a harmonic function w, now satisfying the boundary con- 


ditions 
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“J 
© 2) 


(hi(w, —w), on S,, 
wn; = 90, on S,, 

| , e 

Ao(wo — w), on SS, 


where w, and wy, are constant inlet and outlet ambient temperatures and h, and hy are 
positive, constant heat transfer coefficients. 

These boundary conditions approximate a wide variety of modes of surface heat 
transfer, including radiation, convection, and contact resistance. For large h the bound- 
ary conditions approach the Dirichlet conditions satisfied by « and for small h (with 
some reservations as to the character of the limiting process) the Neumann conditions 
satisfied by v. We therefore expect that in this case the resistance will have a value 
t(u) and R(v), namely 


Riu) < Rw) < Ro), (3.1) 


intermediate between 


where 
R(w) = ((w); — (w)o)/Q(v), 


according to our general definition of resistance. In the next section we obtain upper 
and lower bounds on R(u) and R(v) which we then use to derive the inequality (3.1). 
4. Upper and lower bounds on resistance. The problem of bounding the resistances 
R(u) and R(v) is equivalent to bounding the corresponding Dirichlet integrals. Such 
bounds may be obtained in a variety of ways (see, e.g., [2], [3], [4]). For our purposes, 
the most direct approach is the use of Schwarz’s inequality, as suggested by Diaz and 
Weinstein [1]. We first state the results and the conditions on the bounding functions 
and then give a brief derivation of one of them. 
We have 
(uy — Uo)’ /D(f) < Ru) < Dip) /Q*(p,), (4.1) 


(1 — (o)?/Df) < RO < D(p)/Q, (4.2 


where, as before, 


+ 


Dif) = f.ifidV, Do) =| pp.av, 
JV V 


« 
. 


Op) = | pin; dS, 


+ 


ts [ f dS, (f\o = So' | f dS. 
JS; So 


d 


f)s 


II 


In these inequalities f is any non-constant scalar function, continuous on V + S, piece- 
wise continuously differentiable in V, and satisfying any boundary conditions of Dirichlet 
type given on S. Thus in (4.1) 
_ ju, on 8, 
f= 4 
Up, on So, 
while in (4.2) it is unrestricted. The vector p, is any non-zero, solenoidal vector function, 
having a continuous normal component across any surface in V, piecewise continuous 
derivatives in V, and such that its normal component satisfies any given Neumann 
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conditions on S. Thus in both (4.1) and (4.2) 
p:.1 = 0, in V, 
es 


| 
—) 
g 
—s 


Pn; = 
while, in addition in (5.2), 
(Q/S on S 
/ 1» 1 
pn = 4 ‘ 


\—Q/So, on Sp 


As a sample, let us derive the upper bound on R(u). For two arbitrary vector functions, 
p; and q; , Schwarz’s inequality takes the form 


[D(p; , a)’ < Dp) D(q), (4.3) 


where, using obvious notation, 
D(p, , 9:) = [ pa av. 
JV 


If we set u.; , this becomes 
[D(p, , u,)]’ < Dip) Dw). 
Green’s identity gives 
D(p; ,u,;) = [ upn; dS = (uy — Up) Q(p,), 
using the boundary conditions on p; and u. Thus 


(u, — U)*/D&u) < Dp,)/Q*(p,), 


which is the required inequality. The remaining bounds follow in so similar a fashion 
that there is no need to reproduce their derivation here. 

Now note that we may set p; = w,, in the upper bound on R(u) and f = w in the 
lower bound on R(v), giving 


Riu) < Rw) Dw)/Qw)((w), — (w)o), 
Rw) => Rw)Q(w)((w), — (w)o)/D(w). 
We now show that 
D(w) < Q(w)((w), — (wo); 
from which the inequality 
Ru) < RW) < RO (4.4) 
follows. 


Green’s identity gives 


Q(w) = [ w.n; dS = h,S,(w, — (w),) 
8S; 


= -| wn; dS = —hSp(wo — {w)o), 
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and 


’ 


D(w) = h,S,w,(w), + AoSowolw)o — hy | w dS — ho | w dS. 
Se 


/#S, 


Schwarz’s inequality gives 


S, | w dS > || w as| 
8S, JS, 


hi [ w? dS > h,S,(w)? 
8; 


and similarly for the integral over Sy , so that 
Dw) < hi S(w)(wi — (w)1) + hoSo(w)o(wo — (w)o) = Q(w)((w), — (wo), 


which is the desired inequality. 

5. Conclusion. We have now shown that the thermal resistance, i.e., the average 
inlet-outlet temperature difference for unit input heat power, in the case of the Robin- 
Neumann problem is bounded below by the resistance R(u) for constant inlet and 
outlet temperatures and above by the resistance R(v) for constant inlet and outlet 
heat flux. We have also shown how upper and lower bounds on R(u) and R(v) may 
be obtained. Thus we have a method for estimating the temperature rise which is rela- 
tively independent of the specific nature of the boundary conditions. It is felt that this 
method of obtaining broad bounds on temperature rise in many cases is much more 
practical than, for example, the calculation of an exact solution. It is particularly suited 
to the kind of problem encountered, for example, in preliminary design work, where 
one knows little about the specific nature of the heat flow and where one may want 
to evaluate the effect of changes in many design parameters. 

From the mathematical point of view, the foregoing results are quite trivial. It is 
interesting to note, however, that they provide another illustration of the complementary 
nature of the Dirichlet and Neumann conditions. According to the Dirichlet principle, 
among all functions having given constant values on inlet and outlet, the harmonic 
function u makes the Dirichlet integral D(u) smallest. If we relax this condition and 
require only that the admissible functions have a given difference of average values 
between inlet and outlet, the Dirichlet integral is minimized by the harmonic function v. 
The first result gives the lower bound on R(u), the second the lower bound on R(v). 
The upper bounds may be similarly described, using the Thomson principle concerning 


solenoidal vectors. 
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